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ABSTRACT 


An experimental investigation of the transition process on flat-plate 
and concave curved-wall boundary layers for various free-stream turbulence 
levels was performed. Where possible, sampling according to the 
intermittency function was made. Such sampling allowed segregation of the 
signal into two types of behavior— laminar-like and turbulent-like. Results 
show that for transition on a flat-plate, the two forms of boundary layer 
behavior, identified as laminar-like and turbulent-like, cannot be thought of 
as separate Blasius and fully-turbulent profiles, respectively. Thus, simple 
transition models in which the desired quantify is assumed to be an average, 
weighted on intermittency, of the theoretical laminar and fully turbulent 
values is not expected to be successful. Deviation of the flow identified as 
laminar-like from theoretica 1 laminar behavior is shown to be due to 
recovery after the passage of a turbulent spot, while deviation of the flow 
identified as turbulent-like from the fully-turbulent values is thought to be 
due to incomplete establishment of the fully-turbulent power spectral 
distribution. Turbulent Prandtl numbers for the transitional flow, computed 
from measured shear stress, turbulent heat flux and mean velocity and 
temperature profiles, were less than unity. For the curved-wall case with low 
free-stream turbulence intensity, the existence of Gortler vortices on the 
concave wall within both laminar and turbulent flows was established using 
liquid crystal visualization and spanwise velocity and temperature traverses. 
Transition was found to occur via a vortex breakdown mode. The vortex 
wavelength was quite irregular in both the laminar and turbulent flows, but 
the vortices were stable in time and space. The upwash was found to be more 
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unstable, with higher levels of u' and u'v', and lower skin friction 
coefficients and shape factors. Turbulent Prandtl numbers, measured using a 
triple-wire probe, were found to be near unity for all post-transitional profiles, 
indicating no gross violation of Reynolds analogy. No evidence of 
streamwise vortices was seen in the high turbulence intensity case. It is not 
known whether this is due to the high eddy viscosity over the entire flow 
which reduces the turbulent Gortler number to stable values and causes the 
vortices to disappear, or whether it is due to an unstable vortex structure. 
Predictions based on two-dimensional modelling of the flow over a concave 
wall with high free-stream turbulence levels, as on the pressure surface of a 
turbine blade, would seem to be adequate. High levels of free-stream 
turbulence superimposed on a free-stream velocity gradient (which occurs 
within curved channels) was found to cause a cross-stream transport of 
momentum within the "potential core" of the flow. The total pressure 
within the "potential core" can thus rise to levels higher than that which 
occurs at the inlet to the test section. 

Documentation is presented in two volumes. Volume I contains the 
text of the report including figures and supporting appendices. Volume II 
contains data reduction program listings and tabulated data. 
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NOMENCLATURE 


A Area 

b Bar width 

Cf Skin friction coefficient 

Cp Specific heat or static pressure coefficient, depending on context 
G Gortler Number (p. 4) 

Gr Grashof number 

H Shape factor (=81/82) 

k Conductivity or an empirical constant, depending on context 
i Current 

N Gortler number based on boundary layer thickness 
P Production of shear stress 

Pr t Turbulent Prandtl number 

PSD Power spectral density 

4 Heat flux per unit time and area 

R Wall radius of curvature or resistance, depending on context 
r Local radius of curvature 

Re Reynolds number 

St Stanton number 

T Temperature, mean value 

t Time or instantaneous temperature, depending on context 
TI Turbulence intensity (u'/Uoo x 100) 1 

U Mean stream wise velocity 

u Instantaneous streamwise velocity 

u't' Streamwise turbulent transport of heat, time averaged 
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u ' v ’ Turbulent shear stress, time averaged 
V Voltage 

v't' Cross-stream turbulent transport. of heat, time averaged 
x streamwise distance 

y cross-stream distance 

z cross-span distance 

8 Boundary layer thickness based on 99.5% of freestream velocity 

81 Displacement thickness (p. 16) 

82 Momentum thickness (p. 16) 

e Emissivity or ratio of resistances, depending on context 
y Intermittency (p. 25) ... 

A Non-dimensional wavelength (p. 104) 

X Wavelength of vortices 

v Kinematic viscosity 

II Coles wake parameter 

p Density or autocorrelation (p. 130), depending on context 

t Shear stress or time delay, depending on context 

Subscripts 

c critical or computed value, depending on context 

e eddy value 

h heater 

1 laminar 

p local potential flow value 

pw potential value at wall 

t turbulent 
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tr transition 

w wall value or wall, depending on context 

x based on streamwise distance 

x shear value 

<*> free-stream value 

Superscripts 

' fluctuating component, rms 
+ wall coordinates 

— (overbar) mean value 
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CHAPTER 1 
Introduction 


Despite the attention of many investigators, understanding of the 
boundary layer transition process remains elusive. The sensitivity of 
transition to many factors (free-stream acceleration, the level of free-stream 
turbulence and its characteristics, surface roughness, surface curvature, 
surface heating, wall suction, compressibility and unsteadiness, to name a 
few) makes prediction of the transition process in machines such as gas 
turbines very difficult. Although a very few instances of direct solution of the 
Navier-Stokes equation have recently emerged that yield an extremely 
valuable window to the physics of bypass transition, much of the insight into 
this complex phenomenon is gained by experimentation. The transition 
process is sufficiently complex that observations muse first be made in simple 
geometries with few effects. Later, as understanding builds, more effects can 
be added and more realistic geometries can be investigated. 

1.1). Flat-plate Transition 

The purpose of the first portion of the experimental program was to 

document the effects of three levels of free-stream turbulence on flat-plate 

transitional boundary layers. The experiments provide support for the testing • 

and development of transition prediction models. Specifically, the 

applicability of intermittency-based transition models first proposed by 

Dhawan and Narasimha (1958) which assume a Blasius-type flow for the 

laminar portion and a fully turbulent flow for the turbulent portion were 

tested. Although quite a few researchers have studied the flat-plate transition 

process (see Wang— 1984 for a good review), only a few have used conditional 
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sampling on intermittency to look at the laminar and turbulent portions of 
the transitional boundary layer separately. No previous researchers to the 
authors' knowledge have directly measured the turbulent heat flux and 
Prandtl number in the turbulent part of the intermittent boundary layer. 

Transition on flat-plate boundary layers on smooth walls occurs via a 
T-S (Tollmien-Schlichting) path (Figure 1.1) or a bypass mode depending on 
the free-stream turbulence level (Morkovin—1977). For low free-stream 
turbulence levels, the instability is first manifested in the formation of two- 
dimensional Tollmien-Schlichting waves, which then become unstable in 
the spanwise direction and form hair-pin vortices (Klebanoff, Tidstrom and 
Sargent— 1962, and Perry, Lim, and Teh— 1981). Breakdown to turbulence 
occurs shortly afterwards through the formation of turbulent spots, first 
discovered by Emmons (1951). The bypass mode of transition occurs at higher 
free-stream turbulence levels, turbulent spots forming without T-S wave 
amplification. Spot formation is characterized by a sudden explosion of the 
boundary layer from laminar to turbulent flow (Suder, O’Brien and 
Reshotko— 1988). 

Dhawan and Narasimha (1958) were the first to view the transition 
process as being composed of a Blasius profile alternating with a fully 
turbulent log-law profile. They proposed calculating boundary layer 
parameters within transition by weighting, on the intermittency, the 
corresponding parameters in the fully laminar and fully turbulent flows, each 
at it’s appropriate Reynolds number. Their measurements supported their 
hypothesis. A series of experiments by other researchers in which conditional 
sampling techniques were used to measure quantities within turbulent spots, 
however, generated conflicting results. Wygnanski, Sokolov and Friedman 
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(1976) and Blair (1991) found that the turbulent-zone mean velocity profiles 
along the centerline of the turbulent spot agreed with the log-law, indicating 
that the modelling of Dhawan and Narasimha (1958) could be used for 
transition. Antonia, Chambers, Sokolov, and van Atta (1981), however, 
found that the skin friction required to make the velocity profiles agree with 
the log-law were unreasonably high. Similar results were found by Cantwell, 
Coles and Dimotakis (1978). Mauter and van Atta (1986) found that the shear 
stresses in the center of turbulent spots were 10% to 15% lower than those 
found by Wygnanski, Sokolov, and Friedman (1976), and concluded that 
turbulent correlations can be used for qualitative, but not quantitative, 
descriptions within the spot. A flow visualization study by Gad-El-Hak, 
Blackwelder and Riley (1981) found that the flow in the forward overhang of 
the turbulent spot was relatively passive, being cut off from the bursting 
mechanism at the wall. 

The effects of elevated free-stream turbulence on transition have been 
studied by van Driest and Blumer (1963), Hall and Gibbings (1972), Abu- 
Ghannam and Shaw (1980), Blair (1982), and Wang, Simon and 
Buddhavarapu (1985). The results of a large number of transition 
experiments were examined by McDonald and Fish (1973), who formulated a 
quantitative model which allowed prediction of the onset and extent of 
transition as a function of free-stream turbulence. Blair (1982) and Wang, 
Simon and Buddhavarapu (1985), who both measured the heat transfer in 
transitional boundary layers, found that the temperature profiles lagged the 
velocity profiles and that the turbulent Prandtl number was somewhat 
greater than unity, as deduced from mean profile measurements. 
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1.2). Transition on Concave Walls 


In the second portion of the experimental program, the effects of 
concave curvature on transition were documented. The Taylor-Gortler 
vortices which form on the concave wall (first predicted by Gortler— 1940, see 
Figure 1.2) hasten the transition process by producing unstable cross-span and 
cross-stream inflection point velocity profiles. 

Clauser and Clauser (1937) and Liepmann (1943) were the first to look 
at curvature effects on transition. Both researchers concluded that concave 
curvature had a destabilizing effect on the flow, transition occurring earlier 
than on a flat plate. Wortmann (1969), in a flow visualization study, 
identified three modes of instability. The formation of Gortler vortices was 
the primary instability. The secondary instability manifested itself as a tilting 
of the vortex structure, resulting in highly unstable double inflection point 
velocity profiles. A third order instability in which the vortex structure 
oscillated was then observed. Bippes (1978) also observed a meandering of the 
vortex structure prior to breakdown to turbulence. The formation of vortices 
was found to be described by the Gortler number, G, given by 



The critical Gortler number (Gc) was found to range from 6 to 10, in 
agreement with other researchers, with Gc decreasing with increasing free- 
stream turbulence intensity. Pressure gradients in the direction of the flow 
had little effect on the stability. Swearingen (1985), using smoke visualization 
and hot-wire rakes, found that the breakdown of vortices occurs via two 
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modes— a horseshoe vortex mode and a sinuous mode. Breakdown to 
turbulence, which eventually destroyed the coherent three-dimensional 
structure of the vortex held, occurred shortly afterwards. Inflection points in 
the spanwise direction were found to be more unstable than inflection points 
in the cross-stream direction. McCormack, Welker and. Kelleher (1970), who 
studied the effect of Gortler vortices on heat transfer in a duct, found Nusselt 
numbers 30% to 190% greater on the curved wall than the corresponding flat- 
plate values. This result disagreed with their conclusion, from theoretical 
linear stability calculations, that there should be no net heat transferred due 
to the vortices. 

The effects of concave curvature on turbulent boundary layers is well 
documented. One of the first to study this was Tani (1962), who proposed 
replacing the molecular diffusivity in the Gortler number with the eddy 
diffusivity to obtain a turbulent Gortler number. So and Mellor (1975) found 
a system of longitudinal vortices that were unstable, and that broke up 
downstream, resulting in high turbulence levels. The mean flow was not 
homogeneous in the cross-span direction. Ramaprian and Shivaprasad (1977) 
found the outer region of the boundary layer to be very sensitive to wall 
curvature. The outer region reached a self-preserving form very soon after 
entry into the curve. Mean profiles were found to agree with the log-law, the 
extent of the turbulent core being increased by concave curvature. Shizawa 
and Honami (1983) found similar results. Coles profile parameter (n) was 
found to decrease to zero and even become negative. In a later paper 
(Shizawa and Honami-1985), they suggested that the Gortler numbers may be 
reduced to the stable regime if the eddy viscosity becomes large enough, 
causing any vortex structure within the boundary layer to disappear. Barlow 
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and Johnston (1988 a,b) found, using LIF flow visualization, longitudinal 
vortices that appeared and disappeared randomly in space and time. When 
vortex generators were placed upstream of the curve, however, the vortex 
motion stabilized. Inflows were found to suppress the bursting process, 
outflows to enhance it. Although lower velocities near the wall at the 
outflows would suggest a decrease in the local skin friction, the increase in 
bursting seemed to compensate for the drop, resulting in a relatively constant 
skin friction across the span. They felt that a two-dimensional simulation of 
the flow would be sufficient. Similar conclusions were reached by Simonich 
and Moffatt (1982) in a heat transfer study in which they found that the 
Stanton number varied by only 15%, even under the most energetic inflows. 

1.3). The Measurement Program 

In the tests, a boundary layer is allowed to undergo transition naturally, 
becoming a fully turbulent boundary layer by the end of the test section. The 
effects of three levels of free-stream turbulence were investigated (nominally 
0.32%, 1.79%, and 8%). The wall curvature used in the curved wall cases was 
R=0.97 m. The measurements consist of the following quantities: 

1) . Mean and fluctuating components of streamwise velocity. Mean 

and fluctuating velocities were measured using a horizontal hot-wire 
(TSI Model 1218 Boundary Layer Probe) in isothermal flows. Free- 
stream fluctuating components were measured using a special rotating 
slant wire. 

2) . Mean temperature profiles. A thermocouple probe consisting of 

butt-welded 76 pm (3 mil) dia. chromel-constantan wires held between 
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two hypodermic needles, as described by Wang and Simon (1987), was 
used. 

3) . Local Stanton number. Thermocouples were embedded in the 

test wall for this purpose. 

4) . Shear stress profiles and profiles of the fluctuating component of 

cross-stream velocity were documented in isothermal flows where the 
boundary layer was sufficiently thick. A cross-wire probe (TSI Model 
1243 Boundary Layer Probe) was used for these measurements. 

5) . Intermittency. A horizontal hot-wire was used to determine 

whether the flow was laminar-like or turbulent-like. An analog 
intermittency function was generated so that processing based upon the 
state of the flow (laminar-like or turbulent-like) could proceed as 
appropriate. 

6) . Profiles of the turbulent heat flux v't' were made where the 

boundary layer was sufficiently thick. A triple-wire probe developed 
for this purpose is described below. Measurements of the turbulent 
Prandtl number were made using this probe. 
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Fig. 1.1 Schematic of the transition process on a flat-plate boundary layer at 
low turbulence intensities. (From Schlichting-1979) 



Fig. 1.2— Schematic diagram of Gortler vortices. 
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CHAPTER 2 


Experimental Apparatus, Techniques, and Qualification 
2.1). Experimental Apparatus 

Wind Tunnel . A schematic of the test facility is shown on Fig. 2.1. 
The tunnel was originally designed and built by Wang (1984). It has since 
been modified by adding two settling chambers, a honeycomb pack and a 
motor controller. Capability of increasing the free-stream turbulence to ~9% 
has also been added. A brief description of the facility follows. 

The wind tunnel is of the low-speed, open-return type. Air is drawn 
through filters capable of filtering particles larger than 5 pm which could 
damage the fine hot-wires downstream. A 7 hp centrifugal blower with a 
rated capacity of 5500 cfm and driven by a 3-phase, 230 V, 10 hp motor forces 
air through a series of grids and a honeycomb section. The grids provide 
resistance, aiding in the redistribution of the flow, while the honeycombs 
remove swirl and orient the flow axially. An oblique header deflects the flow 
into a heat exchanger used to control the flow temperature. A 3/4 hp motor- 
driven centrifugal pump circulates water from a 40 gallon tank which serves 
as a thermal capacitor. Immediately downstream of the heat exchanger is a 
honeycomb pack which re-orients the flow axially. Five screens downstream 
of the honeycomb break up the flow, after which it enters a 10.6:1 contraction 
nozzle. The nozzle has an exit flow aspect ratio of 6:1 to minimize secondary 
flow effects in the test section. The velocity at the exit of the nozzle can be 
continuously varied from 6 m/s to 35 m/s using a motor controller (Louis- 
Allis Lancer Jr. VT, 10 hp). Strong suction was applied at the exit of the 
contraction to re-start growth of the boundary layer. A 2 hp centrifugal fan 
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draws flow from the suction box and returns it to the filter box. The suction 
rate is controlled by a valve in the ductwork. 

Two levels of higher free-stream turbulence in the test section can be 
achieved with the use of inserts. A removable grid designed to generate 
-1.5% free-stream turbulence in the test section could be placed between the 
screen pack and the contraction nozzle. The grid consisted of 2.5 cm 
aluminum strips riveted together in a square array on 10 cm centers. Free- 
stream turbulence intensities of -8% at the inlet of the test section could be 
achieved using an insert section after the contraction nozzle. This insert, 
shown on Fig. 2.2, consists of a bi-plane grid of 4.2 cm OD PVC pipes on 10.8 
cm centers and a 96.5 cm long establishment region to allow for turbulence 
development. The grid was similar to that used by O'Brien and vanFossen 
(1985). Mean and fluctuating velocity measurements at the exit of the 
establishment region (just before the test section) shown on Fig. 2.3 measured 
with a hot-wire indicate velocities that are uniform to within 3%, and 
turbulence intensities that are uniform to within 6%. A rotating slant wire 
(see Russ-1989), used to measure all three components of velocity, showed 
that u'~1.06v' and u'~w’. The turbulence was, thus, quite isotropic. 

A schematic of the test wall is shown in Fig. 2.4. The design is similar 
to that of Wang (1984), the main difference being that the stainless steel/3-M 
P-19 film has been replaced by a lexan/liquid crystal sheet. The lexan allows 
the wall to be bent into a concave configuration without the waviness that 
would have resulted with a stainless steel sheet. Measurement of the 
emissivity of the liquid crystal sheet eliminated the need for the P-19 
reflective film that had been used by Wang (1984) for radiation control. The 
thermal conductivity of the lexan/liquid crystal was also measured so that the 
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temperature drop across this element could be calculated. The measurement 
procedures are described in the Appendix. 

The liquid crystal was added so that the transition process could be 
visualized. Cholesteric liquid crystals change color with temperature, easily 
enabling the temperature field to be mapped. Also, with a constant heat flux 
boundary condition, isochromes correspond to lines of constant heat transfer 
coefficient. The transition location can thus be determined by gradually 
increasing the wall heat flux and seeing where the liquid crystal first changes 
color (this is the point of lowest heat transfer coefficient). The liquid crystal is 
also useful for visualizing variations in wall temperature caused by 
longitudinal vortices in the concave configuration. Cooler lines correspond 
to the downwash between vortices, with hotter lines corresponding to the 
upwash. The formation and growth of the vortices can be monitored, and 
their spacing determined. 

2.2). Instrumentation 

Laboratory Computer . A Hewlett Packard Series 200 Model 16 personal 
computer with 1.3 megabytes of memory, a math co-processor, and Basic 
compiler was used for data reduction and as a controller. The computer is 
linked via an IEEE interface with an external dual disk drive (HP9122) which 
utilizes 3.5 in. double sided diskettes, and an inkjet printer (HP Thinkjet). 
High level HP Basic is the programming language used. 

Hot-Wire Anemometer and Probes . A four channel constant- 
temperature anemometer (TSI IFA-100) was used to drive the hot-wires. The 
anemometer features a built-in microprocessor with non-volatile memory 
which monitors and stores set-up parameters for future use. Built-in signal 
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conditioners enable tailoring of the output signal, maximizing digitizing 
resolution of the A/D converters. 

Four types of hot-wire probes were used. Single wire horizontal wire 
probes (TSI 1218 Boundary Layer probe) were used to measure the mean and 
fluctuating components of streamwise velocity in isothermal flows and in 
near-wall situations. The probe prongs are bent at right angles to the probe 
holder such that prong/flow interference is minimal. Two-wire X-type 
probes (TSI 1243 Boundary Layer "X" probe) were used to measure the 
turbulent shear stress, u'v'. The prongs of this probe are also bent at right 
angles to the probe holder. A specially made triple-wire probe, described 
below, was used to measure the turbulent heat flux, v’t'. 

Thermocouple Probe . A thermocouple probe constructed following 
the design of Blackwell and Moffatt (1975) was used to measure temperature 
profiles within the boundary liyer. Details of the design are presented in 
Wang (1984). A short description follows. Chromel-constantan 0.076 mm 
dia. (3 mil) thermocouple wire butt-welded at their junction was held 
between two supports which were separated by 13 mm and made of 22 gauge 
stainless steel hypodermic needles. The supports are electrically insulated 
from one another. The support/ thermocouple assembly can be rotated 
slightly to align the probe parallel to the test wall and perpendicular to the 
flow. A slight bow in the wire allowed the junction to be placed very close to 
the test wall. 

Pressure Transducer . The transducer is of the variable reluctance type 
(Validyne DP45), and is designed for differential measurements of extremely 
low pressure differences. The diaphragm has a pressure range of 0-8.9 cm (0- 
3.5 in. H 2 O) with an accuracy of 0.5% of full scale. The response was found to 
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be very linear throughout the entire operating range. Calibration of the 
transducer was performed against a micro-manometer (Dwyer Microtector) 
with an accuracy of ±0.06 mm H 2 O. The calibration curve was obtained by a 
linear fit to the data. 

Carrier Demodulator . The analog output of the pressure transducer is 
sent to a demodulator (Validyne CD-15) which provides a stable DC output. 
The gain and offset on the demodulator enabled tailoring the signal to 
maximize digitizing resolution. 

A/D Convertor . The A/D converters consisted of an HP 3437A system 
voltmeter, a Fluke 8840A multimeter, and two Norland (now Hi-Techniques) 
Prowler digital oscilloscopes. The first is a 3 1/2 digit successive 
approximation digitizer capable of sampling up to 3600 samples per second 
(ASCH mode). The Fluke is a 16 bit digitizer. The Prowler is a 12 bit digitizer 
capable of simultaneously sampling two channels of data at speeds up to 100 
kHz. Two buffers store up to 4096 data points each. The two Prowlers may 
also be connected as master and slave, enabling four channels of data to be 
taken simultaneously. Simple data processing can be done within the unit 
using the built-in math, calculus and signal processing routines, or the data 
may be sent via the IEEE-488 interface bus to a computer. 

2.3). Measurement Techniques 

Spectra . PSD (Power Spectral Density) distributions of the hot-wire 
output voltage were obtained using software provided by Jensen 
Transformers, Inc. The Norland Prowler was used to digitize and store the 
hot-wire output, then the contents of it's buffer were transferred to the 
computer where FFT (Fast Fourier Transform) processing was performed. 
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Due to storage limitations in the Prowler (4K per channel), the spectrum was 
calculated by averaging the FFT of ten separate waveforms. Furthermore, this 
was performed at three different acquisition speeds (20 kHz acquisition with 
low-pass filtering of the hot-wire signal at 10 kHz, 2 kHz acquisition with 
filtering at 1 kHz, and 200 Hz acquisition with filtering at 100 Hz) to obtain a 
realistic spectrum over a wide range of frequencies. A total of 30 traces was 
therefore required to obtain a PSD. Due to the long times required, only one 
PSD (taken in the free-stream at the test section entrance) was obtained for 
each case. - 

Mean and fluctuating velocity . Mean velocities were obtained by 
averaging the instantaneous velocities measured using a hot-wire oyer 
approximately 30 seconds and 5000 data points. Digitizing was performed 
using the HP 3437A single channel A/D meter. Fluctuating velocities (rms) 
were found according to the formula 



The resultant mean velocity and turbulence intensity profiles were smooth. 

The wall position was found in two steps. First, the probe was placed 
close to the wall, and traversed towards the wall in 50 pm increments until 
the hot-wire output voltage ceased to change. This meant that the prongs of 
the hot-wire probe were in contact with the wall. The probe was then 
traversed away from the wall until an abrupt increase in the anemometer 
output voltage was observed. This position was taken to be the y=0 position. 
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Although this sounds risky, no wire ever broke during this procedure. The 
calibration curves were found to be very repeatable, indicating that no 
straining of the wire took place. The second step in finding the wall position 
was performed during data reduction after the raw data was acquired and 
stored on disc. In turbulent flows, the law-of-the-wall, given by 

u+=2.44 lny+ + 5.0 (2-2) 

u+ _ u / u pw y + _ yupwVQ / 2 
where U >/Cf7p f v 

' / * * 

was used to find the local skin friction coefficient (Cf) as well as the correction 
on the y-position. This is the Clauser (1956) technique. The data was fit to the 
equation u+=y+ in laminar flows, upstream of transition, to find the local skin 
friction. The y=0 position was found by linearly extrapolating the velocity vs. 
position data to the wall. Within transition, no comparable technique is 
available, so skin friction values were determined by fitting the near-wall data 
points to the u+ vs. y+ curve, while limiting the range of the y-corrections to 
those obtained in. the laminar and turbulent regimes (typically 50 pm). The 
skin friction values thus obtained were checked by a momentum balance. 
Agreement was typically within 20%. 

Once the proper y-corrections were made, displacement thickness ( 8 i) 

and momentum thickness ( 82 ) were calculated according to their definitions: 

S| 

/ o u p d y = ; o (u p -u)dy 
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(2.4) 


where 



-u)dy 


Up = u_ 


(2.5) 


for the flat-wall cases and 



u 


pw 


l-y/R 


( 2 . 6 ) 


for the curved-wall cases. Up W was found by fitting equation 2.6 through two 
points in the free-stream. R is taken to be positive for concave curvature. 
Substituting equations (2.5) and (2.6) into (2.3) and (2.4) yields 


8 . = J„( 1 -Tr:) d >’ 


(2.7) 

( 2 . 8 ) 


for the flat-wall cases and 



(2.9) 



( 2 . 10 ) 


for the curved- wall cases, respectively. Once 8i and 82 are found, the shape 
factor (H=8 i/ 82) and momentum thickness Reynolds numbers can then be 
calculated. 
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Momentum balances were performed using the equation 



d5, 

= "dT + ( 28 ! + ®')' 


du 


p* 


dx 


( 2 . 11 ) 


The balance is rough, as the term d82/dx cannot be accurately determined 
due to the sparsity of 62 values. 

The evaluation of uncertainties in hot-wire measurements is very 
difficult. Even if the uncertainty in the hot-wire calibration is made 
arbitrarily small, there is always the doubt whether the hot-wire response 

inferred from a static calibration is applicable over the frequency range of 

/ * * . 

interest. Perry (1982) states that errors as high as 10% in the mean square 
energy distribution of the turbulence are possible, but that the broad-band 
turbulence results are much less affected since the energy containing 
components of the turbulent motions are mainly weighted toward the low- 
frequency end. The reader is referred to Perry (1982) for further discussion of 
uncertainty in hot-wire measurements. The uncertainty of the hot-wire 
measurements in this thesis will be taken to be 5% for the single-correlation 
measurements (u’ and v') and 10% for the cross-correlation measurements 
(u'v' and v't’). These values are consistent with the scatter in the 
measurements as observed by the author. 

Shear stress . The methodology for measuring shear stress is given in 
Buddhavarapu (1984), and will not be repeated here. In contrast to 
Buddhavarapu (1984), however, where the rotating hot-wire technique was 
applied, the present measurements were taken with a cross-wire probe where 
data is available from both sensors simultaneously. It is digitized using the 
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Norland Prowler. Data points (4k) were taken at a 50 Hz sampling rate. 
Procedures for calibration and alignment of the probe with the flow are 
described in Kim (1986). The shear stress data was normalized by both U pw 
and U t (the shear velocity). Fluctuating quantities (u* and v') were 
normalized on U pw only. 

Stanton number . Stanton numbers (St) were measured using 
thermocouples embedded in the wall. The thermocouples are spaced 2.54 cm 
(1 in) apart along the centerline of the test section. Thermocouples are also 
spaced 5.08 cm (2 in.) apart in the spanwise direction at stations 1, 5, and 6 
within the center 30 cm (12 in.) span. At stations 2, 3 , and 4/thermocouples 
are spaced at 2.54 cm (1 in.) intervals to provide greater resolution. 
Additional thermocouples were provided to measure the free-stream 
temperature and the temperature difference across the fiberglass insulation. 
An additional lead to the voltmeter was provided to check the voltmeter 
zero-point. 

The thermocouples were routed to an isothermal box where they were 
soldered to copper wires. Heat shrink tubing insulated the junctions. The 
isothermal box consisted of two sets of aluminum blocks of nominally 20.3 
cm x 30 cm x 2.54 cm dimensions which sandwiched the thermocouples. 
Foam was used to seal the edges. The blocks were then wrapped with 
fiberglass insulation. A paper cover isolated the box from room air 
movements. Two additional thermocouples served as ice-bath references, 
and a third thermocouple measured any difference in temperature between 
the two sets of aluminum blocks. A 150 channel scanner (Fluke 2205A) and 
digital multimeter (Fluke 8840A) were used to acquire the thermocouple 
voltages. 
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The Stanton number is a non-dimensional form of the heat transfer 
coefficient which relates the actual heat transferred to the flow to the 
maximum heat-carrying capacity of the flow. The quantities needed to 
determine St are the wall heat flux, the wall temperature, the free-stream 
velocity and temperature, and fluid properties. The free-stream velocity and 
temperature and the fluid properties are easily obtainable. The wall heat flux 
was computed by measuring the voltage drop across a precision resistor (~2 Q) 
to find the current (i^) and the voltage drop across the heater (V^). The 
power dissipated in the heater per unit area ( Q h ) was computed from 



( 2 . 12 ) 


The power factor was measured to be very close to unity. Corrections were 
made for back heat loss (through the fiberglass insulation), streamwise 
conduction and radiation. The conductivity of the lexan/liquid crystal 
composite (k w ) and the emissivity of the liquid crystal (e) were measured as 
described in Appendix A1 and A2. Wall temperatures were computed by 
measuring the thermocouple voltages and the heat flux. Since the 
thermocouples are located behind the lexan/liquid crystal composite, 
corrections must be made for the temperature drop within the composite. 
This correction was computed from 


q w Ax 

k w (2.13) 


where Ax is the composite thickness and q w is the heat generated in the 
heater minus the back heat loss and the streamwise conduction divided by 
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the heated area. This corrected wall temperature is also the temperature used 
in calculating the radiant heat transfer. 

An uncertainty analysis on the Stanton number, performed with the 
wall in the straight-wall configuration, yielded a maximum uncertainty in St 
of 6%. This number was calculated using the methodology of Kline and 
McClintock (1953) whereby the uncertainty of the input parameters to St (8xi), 
all based upon the 95% confidence level, could be combined to yield the 
uncertainty in St (8St) at the same confidence level. This is the root-sum- 
square method given by 


8St = 


'ast' 


1 L 


v 9xj. 


n2 


8xj 



The partial derivatives in the above equation were evaluated by slightly 
perturbing the values of the input parameters one by one and observing their 
effect on St. The use of a computer program made the calculation very easy. 

Mean Temperature Profiles . Mean temperature profiles were 
measured using the thermocouple probe described earlier. The local wall 
temperature was determined by linearly interpolating between the wall 
temperatures obtained from the thermocouples upstream and downstream 
of the probe location. The y-position correction to the temperature profiles 
was obtained by comparing the near-wall data points to the temperature 
gradient line calculated from the wall heat flux, which was measured 
independently. An example is given on Fig. 2.5. 

Data was reduced to wall coordinates (T + vs. y + ) according to the 
equation 
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(2.14) 



P C p 


The wall shear (t w ) is that deduced from the corresponding velocity profile 
measured in the unheated flow. 

Energy Balance . Energy balances were performed by comparing the 
enthalpy thicknesses obtained by integrating the wall heat flux with that 
obtained from the mean velocity and temperature profile measurements at 
each station. The reader is referred to Appendix E of You (1986) for the 
numerical integration scheme employed. , . . 

Triple-wire Probe for Turbulent Heat Flux Measurements . A probe 
developed for this purpose (Kim and Simon— 1988) is shown on Fig. 2.6. In 
the method used for the probe in this particular reference, two parallel hot- 
wires, operated at different overheat ratios, were used to deduce the 
instantaneous temperature in the flow. The extreme sensitivity of this probe 
to spatial averaging and free-stream temperature variations made operation 
of this probe quite inconvenient, however. The low frequency response of 
the probe (estimated at a few hundred hertz) did not present serious 
measurement problems, but a higher frequency response would have been 
desirable. To solve these problems, it was decided to use the same probe 
geometry, but operate one wire in a constant-current mode as a resistance 
thermometer to measure the instantaneous flow temperature, and to operate 
the other two wires as constant-temperature wires in a standard cross-wire 
configuration to obtain the instantaneous u and v velocity components. The 
disadvantage of this method is that a fine (~1 pm diameter), fragile platinum 
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wire must be used to measure the temperature and sensor survivability can 
pose problems. The advantages, though, are that the measurement of 
temperature is fairly direct, with much less sensitivity to spatial averaging 
effects and no sensitivity to slow ambient temperature variations. 
Furthermore, by using frequency compensation, the frequency response of the 
cold-wire can be greatly improved. 

A scheme similar to that used by Hishida and Nagano (1978) for two- 
wire measurements was chosen as the compensation scheme. In this 
method, the heat transfer coefficient over a cold-wire was estimated from the 
parallel hot-wire signal. It was shown in their paper that 



where V c i = voltage across cold-wire if it had an infinite frequency 

response 

Vl = voltage across cold-wire (measured) 

V 21 = (cold-wire current)x(hot-wire resistance) (constant) 

V 2 = voltage across hot-wire (measured) 

e = ratio of hot-wire and cold-wire resistance at a reference 
temperature (constant) 

k = empirical constant determined from a frequency response 
test 
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It is to be noted that in an isothermal flow, dVi/dt =0, and eqn. (2.15) reduces 
to V c i=Vi as it should. The compensated resistance of the wire (Rc), and, 
thereby, the compensated flow temperature may then be obtained from 



where i is the current through the wire. Once the instantaneous flow 
temperature is known, the hot-wire signals can easily be compensated. The 
constant, k, in eqn. (2.15) is determined from a square-wave test in which the 
probe was placed in a flow and the cold-wire current stepped down from 30 
mA (wire is heated) to 1 mA (wire cools to essentially room temperature). 
The voltage across the wire along with it's derivative is monitored on a 
digital oscilloscope as it cools in a characteristic exponential fashion. A 
compensated voltage variation is calculated from eqn. (2.15) using a guessed 
value of k. The k-value that yields a step drop in V c i is taken to be the correct 
value. The frequency response of the probe, with compensation, was 
measured to be 4 kHz. 

A circuit built for this purpose (Fig. 2.7) consists of a current source, an 
amplifier and a differentiator. A current source of 1 mA drives the cold-wire. 
The voltage across the wire is amplified 200 times then sent through a 
differentiator. The noise requirements on the circuit are tight. The cold-wire 
has a nominal resistance of 50 ohms with the variation of wire-resistance 
within the heated boundary layer being less than 0.05 ohms (rms). For a 
nominal wire-current of 1 mA, this corresponds to only a 50 p V (rms) 
variation. The circuit noise must be much smaller than this to get an 
adequate signal-to-noise ratio. The noise of the amplifier is 0.5 pV rms 
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referenced to input, yielding a signal-to-noise ratio of 100; the signal-to-noise 
ratio at the exit of the differentiator is 30. Careful attention to minimizing the 
potential for ground loops is essential to obtain these values. 

The probe was qualified in a zero pressure gradient flat-plate, two 
dimensional, turbulent boundary layer, with a momentum thickness of 
Ree=1487 and a uniform wall heat flux boundary condition (q w "=l78.9 
W / m2). The boundary layer thickness and free-stream velocity were 0.675 cm 
and 26.6 m/s, respectively. The probe was traversed across the boundary layer 
and measurements of u'v', t', u't' and v't’ were made. Prt values were 
determined from it's definition. 






(2.17) 


Two dual channel digital oscilloscopes (Norland Prowler), wired as 
master and slave, were used to digitize the outputs of the hot-wire bridges 
along with the cold-wire voltage' and it's time derivative. Data was sampled 
at 50 Hz over an 80 second period. An IEEE-488 interface bus was used to 
transfer the contents of the oscilloscope buffers to a computer (HP 9816) for 
storage on disc. Data processing occurred off-line. 

Reduced data were compared with that of Blair and Bennett (1984) and 
Gibson, Verriopoulos and Vlachos (1984). All profiles were in excellent 
agreement with the data of these researchers. Measured Pr t values are shown 
on Fig. 2.8. The scatter in the data, typical of direct Prt measurements, is larger 
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in the outer regions of the boundary layer where transport terms and mean 
gradient terms are weak. 

The uncertainty in measured Prt values can be estimated using the 
methodology of Kline and McClintock (1953) discussed above. The error in 
u'v' is assumed to be 10%. The error in v’t' is also estimated to be 10% based 
on scatter in the data and agreement of near-wall v't' values with measured 
wall heat fluxes. The uncertainties in the gradients of velocity and 
temperature become very large in the outer portion of the boundary layer as 
the gradients become small. For example, for the low TI, flat-wall, station 6 
profile (the data is given in the Appendix), the error in the gradients of 
velocity rise from 12% at y/5=0.395 to 53% at y/8=0.85. Similar errors were 
observed for the temperature gradients. Applying the root-sum-square 
propagation of the uncertainty, the errors in Prt increase from 22% to 85% at 
the above y/8 locations, respectively. . The uncertainty in Prt increases sharply 
for this case, at y/ 8=0.71. The uncertainty of the other Prt profiles is expected 
to behave in a similar manner. 

Intermittencv Circuit . A circuit for determining when the flow is 
laminar or turbulent has been constructed. The output of the circuit, an 
analog signal which is high when the flow is turbulent and low when the 
flow is laminar, is called the intermittency function. The intermittency value 
(y) can be found simply by time-averaging the intermittency function. The. 
intermittency function can also be used to conditionally sample the signals so 
that data is processed only when the flow is laminar-like or turbulent-like. 

The circuit takes advantage of the much larger time derivative of the 
turbulent-like signal as compared to the time derivative of the laminar-like 
signal: the hot-wire-anemometer signal is processed by a series of filters. 
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differentiators and rectifiers (Fig. 2.9). At the level detector, the signal is 
compared to an adjustable threshold value. If it is higher than the threshold, 
the output signal of the level detector is high (turbulent-like). It is low 
(laminar-like) otherwise. The analog signal, thus obtained, can be used to 
conditionally sample other quantities, tagging them to either laminar-like or 
turbulent-like behavior. 

It can be seen in Fig. 2.9 that the intermittency measuring unit has two 
channels— the direct channel and the differentiated channel. The two 
channels are used to solve the problem of zero-crossing. This problem is 
explained in Fig. 2.10 where the time-derivative of a turbulent signal is 
shown. When this signal is compared to a threshold value at the level 
detector, the turbulent flow signal is generally higher than the threshold 
giving a recorded "high" signal. However, it is also seen that the signal 
unavoidably becomes smaller than the thresh Did as it crosses zero even 
though it is known to be from the turbulent-like flow. During this time the 
flow is falsely declared laminar. This is the zero-crossing problem. The 
circuit uses the second derivative of the signal (differentiated channel) to 
correct for this. This differentiated signal retains the characteristics of the 
first-derivative with one important difference. It is high when the zero- 
crossing event takes place. When the second time-derivative is compared 
with the threshold value there will again be regions where the flow is falsely 
declared laminar. .However the times during which each of the two channels 
is at fault generally do not coincide. An "OR" gate is then used to combine 
the two signals. Its output is high when either of the two signals is high and 
is low only when both inputs are low. The number of points falsely declared 
laminar is thus greatly reduced. A high-pass filter at the output of the OR 


26 



gate then eliminates the remaining points falsely declared turbulent. The 
threshold values of the two level detectors are adjustable and are tuned for 
each different flow situation. A tuning procedure that has been found to 
work well has been established. 

An example of the circuit performance in the transition region is given 
on Fig. 2.11. It may be seen that the circuit does a good job of discriminating 
between laminar and turbulent flow. The main deficiency of the circuit is 
that it uses a criterion based on the derivative of a hot-wire signal and not on 
eddy transport (u'v')--the mark of turbulence. The advantage of this 
technique, however, is that the wall can be approached very closely with the 
single-wire probe, something not possible with the bulky cross-wire probes 
necessary to determine u'v'. 

The circuit does not perform as well in discriminating between 
boundary layer and free-stream flow in the boundary layer wake region (also 
an intermittent flow). The turbulent fluctuations in the wake decrease in 
intensity, while local pressure fluctuations cause unsteadiness in the free- 
stream, making a criterion based on velocity difficult to implement. A better 
way of discriminating the two regimes in this flow which are intermittent at 
the edge of the boundary layer is to heat the wall and use a criterion based on 
temperature. This technique assumes that mixing in the turbulent portion of 
the wake is thorough enough such that the flow temperature is higher than 
in the laminar region which is at the free-stream temperature. The 
advantage of this technique is it's insensitivity to velocity fluctuations. The 
disadvantage, of course, is that an additional high frequency response 
resistance thermometer (e.g.- a 1 |im Pt wire) is needed. Also, this technique 
requires heating the wall, which affects the transition start location. 
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A comparison of the velocity and temperature measurement schemes 
is shown on Fig. 2.12, where a hot-wire voltage trace and a signal from a 1 pm 
dia. Pt resistance thermometer are presented. The data was taken using the 
triple-wire probe described earlier. The wake passing seems to be more clearly 
indicated from the cold wire signal than with the hot-wire signal. 

In conclusion, the intermittency circuit is seen to give good results in 
cases where there is a clear distinction between regimes (e.g. - transitional 
boundary layers). The circuit does not perform well in the wake region of 
turbulent boundary layers and, it is suspected, in the outer portion of 
transitional boundary layers as well. 

2.4). Test-wall and Tunnel Qualification 

Mean velocity and streamwise turbulence intensity measurements 
within the potential core of the flow exiting the nozzle showed a peak-to-peak 
variation in velocity of 0.2% about a nominal velocity of 27 m/ s. 
Measurements of mean temperature within the flow exiting the nozzle 
showed a peak-to-peak variation of 0.02°C. The static pressure coefficient (Cp) 
was adjusted to within 1.79% all along the test wall for both low and high TI 
cases. 

Qualification of the test section in a flat wall configuration with regard 
to the transition location was initially performed by heating the wall and 
visualizing transition using a liquid crystal sheet. Transition was assumed to 
occur at the location where the liquid crystal first changes color as the heat 
flux is gradually increased. This corresponds to the highest wall temperature, 
or, since the wall heat flux is essentially uniform, the location of lowest heat 
transfer coefficient. Various parameters such as the leading edge suction flow 
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rate and the suction slot width were optimized such that transition occurred 
as far downstream as possible for a given free-stream velocity and for the low 
free-stream turbulence intensity case. The outer flexible wall was adjusted 
such that there was no pressure gradient along the wall. 

This method of determining the transition start location is not reliable, 
however, due to the destabilizing effect of heating on the boundary layer. It is 
well known (Schlichting-1979) that the heating of a surface in air causes an 
inflection in the near-wall velocity profile due to a local increase in viscosity. 
Transition was expected to occur earlier in a heated boundary layer, and the 
transition length was expected to decrease. The variation of Stanton number 
along the wall (low TI case) for two wall heat 'fluxes taken in the test facility is 
shown on Fig. 2.13. It is seen that while transition occurs over a shorter 
length with increasing wall heat flux, as expected, the transition start location 
is not affected. This is due to encroachment of the side-wall influence 
towards the centerline of the test wall. Transition was observed to occur first 
at the tunnel end-walls, due to possible corner flow effects, then propagate to 
the tunnel centerline. The transition start location was thus fixed by the end 
wall effects, obscuring the influence of heating. Transition start was defined 
in this study as the location where the near-wall intermittency as measured 
in the unheated flow reached 5%. 

The St values measured in the laminar flow (before the onset of 
transition in Fig. 2.13 do show, however, how accurately St can be measured. 
Except for a small unheated starting length effect (the first five points), the 
Stanton numbers in the laminar boundary layer are seen to be in excellent 
agreement with the accepted correlation for a constant wall heat flux 
boundary condition. The slight dip in Stanton number values below the 
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laminar correlation is caused by a decrease in the wall heat flux with 
streamwise distance due to increasing radiant heat loss. A STAN5 (Crawford 
and Kays— 1976) simulation with the measured wall heat flux input as the 
wall boundary condition yielded Stanton numbers 4% lower than the 
constant wall heat flux correlation at the start of transition— in excellent 
agreement with the above trends in the data. 

The Reynolds numbers based on displacement and momentum 
thicknesses at the beginning of transition (unheated flow) were measured to 
be 1920 and 737, respectively. The free-stream turbulence intensity, measured 
using a cross-wire (TSI Model 1243 Boundary Layer Probe) rotated into two 
positions to get all three velocity components, was 0.32%. A plot of the 
Reynolds number based on displacement thickness vs. the free-stream 
turbulence intensity for the present study is shown on Fig. 2.14. Transition is 
seen to occur slightly earlier for the present low TI case than for other 
researchers due to the sidewall influence. 

An energy balance was performed by integrating the wall heat flux 
along the centerline of the ,est wall and comparing this with the increase in 
energy carried in the boundary layer flow as calculated from the mean 
velocity and temperature profiles. The closure was within 3%. 

Further qualification of the test section and measurement techniques 
was performed by comparing data measured in the flat wall transitional flow 
with that of other researchers. Measurements of the mean velocity profiles, 
shape factor (H) and intermittency all were consistent with other researchers’ 
results. For the heat transfer data, it was decided to work with the lowest wall 
heat-flux level which would still give reasonable wall-to-freestream 
temperature differences at the end of transition (the location of smallest 
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temperature difference). The nominal heat flux chosen was 175 W/m 2 , 
which yielded a minimum temperature difference of about 4° C for the 
TI=0.32% case. As mentioned earlier, the transition process was significantly 
affected by the heating at this wall heat flux, the near-wall intermittency 
increasing from 45% with no heating to 98% with heating at a selected point 
within the transition zone. This unfortunately means that a precise 
comparison of the heated and unheated data cannot be made for the low TI 
case. For this reason, no heat transfer data will be presented for the low TI 
case, except for measurements of the turbulent heat flux and turbulent 
Prandtl number (Prt). • Transition for the high TI case (TI=1.79%) was not 
affected by heating, the intermittency remaining invariant with the wall heat 
flux. All heat transfer data will therefore be presented for cases where 
11=1.79% and higher. 

A summary of boundary layer parameters for each of the five cases to 
be discussed is presented on Tables 1 through 5. 
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St. x 

Upw 

6 

61 

62 


Re82 

Cf 

(cm) 

(m/s) 

(cm) 

(cm) 

(cm) 

xlE-5 


xlE*3 

1 0.114 

26.15 

0.141 

0.0369 

0.0155 

1.964 

265-7 

165 

2 0.343 

26.26 

0.216 

0.0692 

0.0275 

5906 

4736 

0.96 

3(L) 0.572 

26.17 


0.0951 

0.0383 

9.748 

6525 

0.71 

3(Tu) 

m 


0.1757 

0.0790 

M 

13^8 

2.70 

3(Tr) 

** 

0.324 

0.0950 

0.0390 

m 

664.8 

0.78 

4A(L) 0.600 

26.06 


0.1120 

0.0461 

12.55 

754 

0.79 

4A(Tu) " 

* 


0.1346 

0.0925 

* 

1450 

315 

4A(Tr) " 

m 

0.694 

0.1107 

0.0553 

m 

867.4 

1.40 

4(L) 0.600 

26.09 


0.1014 

0.0437 

1353 

7395 

0.79 

4(Tu) 

M 


0.1231 

0.0851 

•• 

1438 

3.10 

4(Tr) " 

•i 

1.339 

0.1019 

0.0577 

•t 

9759 

1.90 

5CL) 1.029 

2672 


0.0914 

0.0449 

17.76 

776.6 

1.10 

5(Tu) " 

" 


0.1423 

0.1044 

* 

1805 

362 

5(Tr) “ 

H 

1.010 

0.1224 

0.0860 

m 

1487 

377 

6 1.257 

32.64 

1.105 

0.1437 

0.1054 

24.82 

2060 

370 


Table 1 --Summary of boundary layer parameters for Case 1. Values 
conditionally sampled on laminar, turbulent and transitional flow are 
denoted by (L), (Tu), and (Tr), respectively. 
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Table 3~Summary of boundary layer parameters for Case 3. 



St 

X 

(cm) 

Upw 

(m/s) 

6l 

(cm) 

«2 

(cm) 

Rex 

X1E-5 

Re& 2 

Q 

X1E+3 

ATw 

(C) 

Qw" 

(W/m2) 

1 

0.069 

16.53 

0.0541 

0.0213 

0.917 

2192 

223 

— 


2(d) 

2(U) 

0.356 

0.356 

17.24 

17.23 

0.0317 

0.1160 

00164 

0.0531 

3760 

3-757 

172 9 
561 

4.60 

2.10 

% 

— — 

3(d) 

3(u) 

0.610 

0.610 

17.06 

17.11 

0.1407 

0.1623 

00996 

0.1124 

6.369 

6.403 

1044 

1161 

4.60 

415 

395 

423 

147.6 

147.6 

4(d) 

4(u) 

0676 

0.676 

17.14 

17.13 

0.1532 

0.2467 

0.1167 

0.1620 

9244 

9234 

1231 

1917 

520 

420 

366 

420 

146 3 
1463 

5(d) 

5(u) 

1.130 

1.130 

16.76 

16.76 

0.2436 

0.3679 

0.1696 

0.2716 

11.64 

11.65 

1954 

2601 

4.70 

370 

4.15 

4.31 

146.6 

1466 


Table 4-*Summary of boundary layer parameters for Case 4. Upwash and 
downwash are denoted by (u) and (d), respectively. 


St. 

X 

Upw 

R«x 

Cf 

ATw 

Qw" 


(cm) 

(m/s) 

xlE-5 

xlE*3 

(C) 

(W/m2) 

1 

0.069 

17.70 

0965 

6.00 

309 

2169 

2 

0356 

17.70 

3661 

590 

367 

2136 

3 

0.610 

17.70 

6.635 

530 

4.11 

2114 

4 

0.876 

17.70 

9543 

530 

404 

211.6 

5 

1.130 

17.70 

1234 

5.00 

4.21 

210.6 


Table 5-Summary of boundary layer parameters for Case 5. 
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Fig. 2.1 -Schematic of test facility (plan view). 
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Fig. 2.3— Mean and fluctuating velocity measurements at the exit of the 
turbulence establishment section (from Russ-1989). 
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Fig. 2.4-Cross-section of heated test-wall (not to scale). 
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Fig. 2.5— Illustration of the method used to determine the distance of 
thermocouple probe from the heated wall. 
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Fig. 2.6— Schematic of triple- wire probe used to measure turbulent heat fluxes. 
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2.5 pm dia. platinum plated tungsten 

I/d = 200, 1 = 0.5 mm 

separation between wires: 0.35 mm 
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Fig. 2.8~Comparison of Prj values measured' using the probe to the data of 
Blair and Bennett (1984) and Gibson, Verriopoulos and Vlachos (1984). 



Direct channel 



Fig. 2.9-Block diagram of intermittency measuring unit 
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Time [$] 


Fig. 2.11-Intermittency circuit performance in a transitional boundary layer. 
Position in the boundary layer is as noted. 
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Fig. 2.12— Comparison of signals in the wake of a. turbulent boundary layer 
using a), a hot-wire, and b). a cold-wire. 
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Fig. 2.13— Variation in Stanton number vs. Re x for two wall heat fluxes. 
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CHAPTER 3 
Results and Discussion 


3.1). Case 1 — Flat- wall, TI=0.32% 

Free-stream Turbulence Intensity and Spectra . A power spectral 
density (PSD) distribution vs. frequency of the streamwise velocity 
component measured using a horizontal wire (TSI Model 1218 Boundary 
Layer Probe) is shown on Fig. 3.1.1. The power spectrum has a pronounced 
peak at 29 Hz. This peak has been traced (using an accelerometer and a 
vibration analyzer) to -a rocking motion of the centrifugal blower on its 
mounts, resulting in a slight unsteadiness in free-stream velocity. All 
reasonable effort has been applied to minimize this fan motion. This 
frequency is not expected to influence the transition process as the minimum 
critical frequency for amp’ification of disturbances is estimated from linear 
stability theory to be 1600 Hz. The spectrum is seen to be relatively clean 
otherwise. The comparison on Fig. 2.14 supports the conclusion that the 
effect of this rocking motion on transition is minimal. 

Results of measurements of the free-stream turbulence intensity vs. 
streamwise distance using a cross-wire rotated to two positions (TSI Model 
1243 Boundary Layer Probe) are presented on Fig. 3.1.2. The w’ component of 
turbulence was measured at one station only. It is seen that u' is roughly 
twice the value of either v' or w’, with the values remaining constant all 
along the test section. The low-frequency unsteadiness discussed above is 
expected to be the source of the non-isotropy. The free-stream velocity was 
nominally 26.5 m/s. 
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Intermittency Profiles . The intermittency profiles taken at various 
stations are shown on Fig. 3.1.3. The intermittence is seen to decay to zero as 
y/8 approaches 1, as expected. The profiles at stations 4, 5, and 6 show the 
s am p qualitative behavior, namely, a relatively flat value in the region 0.05 < 
y/8 < 0.15 (the turbulent core) followed by a decay to zero for y/8 > 0.15. This 
decay is due to both the entrainment of the free-stream flow into the 
boundary layer (the wake region) and intermittent turbulent spot passing. All 
intermittency values quoted below correspond to the intermittency values in 
the near-wall region where the profiles are flat. 

Mean Velocity Profiles . Profiles of the mean velocity sampled on 
intermittency at stations 3 to 6 are shown on Fig. 3.1.4. The distance away 
from the wall has been normalized on the boundary layer thickness of the 
transitional flow profile. Two characteristics are immediately apparent in all 
the profiles. First, the turbulent boundary layer is thicker than the 
corresponding laminar boundary layer, as expected, due to bursting and 
subsequent turbulent spot formation. Second, the turbulent boundary layer 
profile is flatter than the corresponding laminar profile, resulting in a cross- 
over between the two. The transitional flow profile is, by definition of the 
intermittency, a composite of the turbulent and laminar profiles, and must lie 
between the two. Evolution of the transitional flow profiles is seen to be 
from laminar-like to turbulent-like between stations 3 and 5. 

Plots of the mean velocity sampled on intermittency and normalized on 
wall coordinates at various locations are shown on Fig. 3.1.5. The transitional 
flow profile is seen to evolve from the Blasius profile to the fully turbulent 
log-law profile as noted above. The velocity profiles sampled on 
intermittency, however, do not agree with either the Blasius or log-law 
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profiles in the transition region. The laminar profile increasingly deviates 
from the Blasius profile as transition proceeds while the turbulent profile 
deviates most from the log-law profile early in the transition process. 

A plot of the local skin friction (Cf) values deduced from the near-wall 
velocity gradient (in the laminar-flow case) or by fitting the near-wall data to 
the log-law (in the turbulent-flow case) is shown on Fig. 3.1.6. The skin 
friction corresponding to laminar flow increasingly deviates from the 
laminar correlation as transition proceeds. The higher stress at the wall is 
believed to be due to disturbances in the laminar flow regime as a result of 
nearby turbulent spot passage. A near-wall hot-wire voltage trace in the 
intermittent region shown on Fig. 3.1.7 illustrates this. Although the 
transition from laminar to turbulent flow is quite sharp at the leading 
interface, the laminar flow requires some time to relax back to a nominally 
laminar state. If the inlermittency is high enough, i.e., spots pass frequently, 
the laminar boundary layer is continually disturbed, resulting in higher 
velocities near the wall (than if there were no disturbance) and, consequently, 
higher Cf values. Values of Cf in the turbulent flow, but at the beginning of 
transition (Fig. 3.1.6), are seen to be lower than the fully-turbulent correlation 
values. This could be due to a less than complete establishment of the full 
turbulence spectrum, i.e., only relatively large eddies are present at this stage 
of the transition process and turbulence cascading and dissipation is not fully 
established. This has yet to be confirmed, however. 

A similar variation is seen for the shape factor (H) as shown on Fig. 3.1.8. 
As transition proceeds, H for the laminar boundary layer increasingly 
deviates from the laminar value of 2.6, indicating an increasingly non-Blasius 
type profile. Similarly, early in transition, H deviates substantially from the 
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high-Reynolds-number nominal turbulent value of 1.4. This further 
illustrates that the laminar and turbulent regimes in the transitional flow 
cannot be thought of as being composites of Blasius and mature turbulent 
profiles. 

Velocity fluctuation . The rms of the streamwise velocity fluctuation 
(stream wise turbulence intensity) at stations 3 to 5 are seen in Fig. 3.1.9. The 
rms of the laminar profile at first increases with axial distance but then 
reaches a peak value of 8% at station 4, flattening out thereafter. The peaks in 
the profiles are seen to occur at roughly 30% of the laminar boundary layer 
thickness for all stations. The peak rms of the turbulent profile is initially 
high (16% at station 3), indicating a high production of turbulence, but then 
decays to a peak value of 8% as, it is assumed, the dissipation in the boundary 
layer increases. Equilibrium is reached by station 6. The transitional flow 
profile exhibits quite unexpected behavior. The profile initially follows the 
laminar profile due to the low intermittency (approx. 5% at station 3), but 
then jumps to a peak value of 17.5% at station 4, a value larger than the peak 
in the corresponding turbulent profile. Much of this behavior is due to 
intermittent "switching" of the flow between the laminar and turbulent 
regimes as turbulent spots pass the probe. This was first shown by Arnal, 
Juillen and Michel (1978). The accompanying change in the mean velocity is 
illustrated in the hot-wire voltage trace of Fig. 3.1.7. The differences in the 
mean velocities in the laminar and turbulent regimes give rise to a rms 
velocity fluctuation which is greater than that of either the laminar or 
turbulent regime. In fact, the deviation of the transitional flow profile from 
an average profile of u’ 2 is given by 
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(3.1) 


u'c 2 = Yu? + (1 - Y)u ', 2 +7(1 - y)(u, - u/ 

The level of turbulence as indicated by the transitional flow profile is thus not 
a good measure of the true turbulent activity in the transitional boundary 
layer. 

Shear stress profiles . The variation in the shear stress u'v' through 
transition is shown on Fig. 3.1.10. The laminar contribution to the shear 
stress throughout the boundary layer is seen to be quite small for all stations 
except station 5 (where the number of samples is small and where cross- 
contamination between laminar and turbulent regimes is significant). A peak 
in each profile is seen to move progressively toward the wall as transition 
proceeds. The fully-turbulent profile is reached by station 6. Although the 
transitional flow profile is between the turbulent and laminar profiles for all 
stations, it also is affected by the intermittent "switching" from lami.aar to 
turbulent flow. The transitional flow profile is, therefore, also not indicative 
of the true turbulent shear stress in the boundary layer. 

Turbulent Heat Flux Measurements . Measurements of the turbulent 
heat flux normalized on the wall heat flux and sampled on intermittency are 
shown on Fig. 3.1.11. This normalization is more appropriate than a 
normalization based on the freestream velocity and wall-to-freestream 
temperature difference due to the inability of assigning an appropriate 
temperature difference during transition for the uniform wall heat flux 
boundary condition. In transition, it is presumed that the wall temperature 
fluctuates as the flow regime switches over a given spot on the wall. A 
potential advantage of the present normalization is that the turbulent heat 
flux at the wall should vary directly with the intermittency if v't’ sampled on 
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the laminar flow is small. Due to the destabilizing effect of heating on 
transition, the intermittency for these profiles corresponding to stations 3 and 
4 are different than for the unheated data presented above. This heat flux 
data is presented not as a comparison to the hydrodynamic data, but because it 
provides insight into the transition process. 

It is seen in Fig. 3.1.11 that a large increase in the turbulent heat flux 
above the wall heat flux occurs within the turbulent spot at station 3. This 
can occur if the cross-stream eddy diffusion of heat increases in the 
streamwise direction at the expense of convection of heat. The triple-wire 
measurements bear this out. The streamwise heat flux ut was found to 
decrease almost an order of magnitude between St. 3 and 4 in the near-wall 
region and remain relatively constant thereafter. Whether the 
measurements sampled on turbulent-like flow for station 3 drops to unity in 
the very near-wall region (nearer the wall than can be measured) is not 
known. It is possible that the wall transfers more energy to the flow during 
the passage of a turbulent spot due to the higher heat transfer coefficient, than 
during the times the flow is laminar-like when the heat transfer coefficient is 
lower. If so, the wall heat flux would be varying with time according to the 
local flow regime. The time-average energy transferred must, of course, equal 
the time-averaged wall heat flux. A positive slope in v't also suggests 
intense mixing of the flow away from the wall. It is also seen that v t in the 
laminar portion of. the transitional boundary layer is not zero. This does not 
mean that a turbulent transport of heat is present in the laminar boundary 
layer, but simply that v' and t' are correlated due to the unsteadiness of the 
flow. Because v't' in the laminar regime is not small, the value of v’t’ at the 
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wall for the transitional flow profile unfortunately does not go to the near- 
wall intermittency, as was anticipated. 

The turbulent Prandtl number (Prt) sampled on intermittency through 
transition is shown on Fig. 3.1.12. The uncertainty in this data was estimated 
at 20%. The data at stations 4 and 6 show Prt values consistently close to unity 
in the inner half of the boundary layer, as would be expected of fully 
turbulent boundary layers. The data in the outer half of the boundary layer is 
not expected to be reliable due to the very shallow gradients of velocity and 
temperature. The data for station 3, however, show a drop in Prt values 
sampled on turbulent flow substantially below unity in the near-wall region, 
implying that the eddy diffusivity of heat increases faster than the eddy 
diffusivity of momentum as one traverses from the wall outward. This is 
contrary to the expectation that heat acts simply as a passive scalar and also 
contradicts the conclusions of other researchers (e.g.— Kuan and Wang-1990, 
Blair— 1991, and Suder, O’Brien and Reshotko-1988). These are the first 
measurements in which Prt was directly measured, however, previous 
conclusions having been inferred from mean profile data. 
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Fig. 3.1.1— Power spectral density distribution of the streamwise velocity 
component (Case 1). 
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Fig. 3.1.2-Measurement of the free-stream turbulence intensity vs. 
stream wise distance (Case 1). 
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Fig 3.1 .3-Intermittency profiles along the wall (Case 1). □ 
A - station 4a, • - station 4, ■ - station 5, 0 - station 6. 
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Fig. 3.1.5-Profiles of mean velocity sampled on intermittency and normalized 
on wall coordinates (Case 1). A * laminar, # • turbulent, 4* - transition. 
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Fig. 3.1.6— Skin friction sampled on intermittency (Case 1). A - laminar, 
• - turbiilent,-f- • transition. 5 j 
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Fie. 3.1.7— Near-wall hot-wire voltage trace in transition illustrating 
different mean velocities between the two regimes and the relaxation of the 
boundary layer after turbulent spot passage (as in "A"). 


62 







H 


Fig. 3.1.8-Shape factors sampled on intermittency (Case 1). A - laminar, 
# - turbulent, + - transition. 
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Fig. 3.1.9— Turbulenoe intensity profiles sampled on intermitten cy (Case 1). 
▲ - laminar, # - turbulent, B - transition. 


64 











y 

5 995 


Fig. 3.1.11-Cross-stream turbulent heat flux sampled on intennittency (Case 
1). Station 3 is within transition (y=12%), stations 4 and 6 are fully turbulent. 
A - laminar, • - turbulent, 4 * - transition. 
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Fig. 3.1.12-Turbulent Prandtl numbers sampled on intermittency (Case 1). 
Station 3 is within transition (y=12%), stations 4 and 6 are fully turbulent. 
A - laminar, # - turbulent, -H- transition. 
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3.2). Case 2 — Flat-wall, TI=1.5% 

Free-stream Turbulence Intensity and Spectra . The Power Spectral 
Density (PSD) for this case is shown on Fig. 3.2.1. A peak, corresponding to 27 
Hz, is again visible. As in the base case, this peak is caused by rocking of the 

fan. The frequency is slightly lower in this case, however, since the fan speed 
was lowered to move transition a reasonable distance from the leading edge. 
Another much smaller peak is seen at 3800 Hz. The source of this peak 
(electronic noise) was determined by computing the PSD with the hot-wire 
mounted in the tunnel, but in a still flow. The two dominant peaks in this 
PSD corresponded to 3800 Hz and 60 Hz (line noise). The free-stream velocity 
. was nominally 16.7 m/s. 

Stanton number . The Stanton number variation through transition is 
shown on Fig. 3.2.2. The first five points are seen (as in the lower TI case) to 
be higher than the laminar correlation due to the unheated starting length 
effect. The two data points just before transition agree with the correlation. 
Increasing the free-stream turbulence is seen to have a strong effect on 
transition onset, transition moving to Re x =3xl0 5 , or about one third the 
value for the base case. A comparison of the onset location with other 
researchers was shown on Fig. 2.14. The agreement in this case is very good. 
The Stanton number variation through transition is consistent with the data 
i of Blair (1982). 

Intermittencv profiles . Intermittency profiles are shown on Fig. 3.2.3. 
Like the profiles for the base case, the intermittency monotonically increases 
as transition proceeds. 

Mean velocity profiles . Profiles of the mean velocity sampled on 
intermittency are shown on Fig. 3.2.4. The laminar profile is seen to deviate 
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quite strongly from the Blasius profile throughout the transition process 
(much more than in Case 1), indicating a large perturbation due to increased 
free-stream turbulence. The turbulent profile, in contrast, agrees with the log- 
law profile from very early in transition. The above trends are reflected in 
the skin friction, Cf, values plotted on Fig. 3.2.5. The laminar Cf values 
deviate strongly from the laminar correlation while the turbulent Cf values 
remain relatively unaltered. There is no drop below the turbulent correlation 
in the turbulent Cf values as was seen in the lower TI case. Similar trends 
are seen in the shape factor (H), Fig. 3.2.6. 

Velocity fluctuation . The rms of the velocity fluctuations is shown on 
Fig. 3.2.7. The most striking feature of these profiles in comparison with 
those of the base case is the large increase in laminar unsteadiness, which 
even exceeds the turbulent profile rms values in some parts of the 
transitional boundary layer (station 3). The high values are consistent with 
the observed trends in Cf and H for the laminar regime. The turbulent 
profiles evolve as in the base case. The peak values of the turbulence 
intensity drop more or less monotonically with increasing intermittency for 
the two cases. A fully turbulent profile is established by station 5. 

Temperature profile . The mean temperature profiles through 
transition plotted in T+ vs. y+ coordinates are shown on Fig. 3.2.8. A smooth 
variation from a laminar-like to a turbulent-like profile is seen, as was seen 
for the mean velocity. The temperature profiles lag the velocity profiles, 
however, as may be seen by comparing the two profiles at station 3 (see Fig. 
3.2.4 and Fig. 3.2.8). The temperature profile is still evolving when the 
velocity profile has assumed a nearly log-law shape. This is consistent with 
the observations of Blair (1982). 
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Shear Stress Profiles . The variation in shear stress u’v' sampled on 
intermittency through transition is shown on Fig. 3.2.9. As for the low TI 
case, the laminar contribution to the shear stress is seen to be small 
everywhere except in the very near wall region. The peak in the turbulent 
flow profiles decreases in amplitude while moving towards the wall as 
transition proceeds. An equilibrium profile is reached by station 5. 

Turbulent Heat Flux Measurements . Profiles of the turbulent heat flux 
are presented on Fig. 3.2.10. As in the lower TI case, a strong increase in the 
turbulent heat flux above the wall heat flux is seen. This peak is seen to decay 
rapidly, the profile achieving what one would expect of a fully turbulent 
profile by station 5. 

Turbulent Prandtl numbers deduced from the measurements are 
presented on Fig. 3.2.11. The values are seen to all be in the vicinity of unity 
for the fully turbulent profiles (stations 3 to 5), while dipping below unity in 
the transitional flow case (station 2), as was seen previously for the low TI 
case. 
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Fig. 3.2.2-Stanton number variation through transition (Case 2). 
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Fig. 3.2.3-Intermittency profile variation through transition (Case 2). 
O * station 2 , A- station 3a,#- station 3, # - station 4, O" station 5. 
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Fig. 3.2.5— Skin friction sampled on intermix 
• - turbulent, -f- - transition. 
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Fig. 3.2.7-Streamwise turbulence intensity sampled on intermittency (Case 2). 
A * laminar, turbulent, + - transition. 
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Fig. 3.2.8— Mean temperature profiles through transition normalized on wall 
coordinates (Case 2). A - station t , O * station 2, A - station3a,# - station 3, 

■ - station 4,0 - stations. Station 3a refers to data taken at station 3, but at a 
reduced free-stream veloctiy. 





Fig. 3.2.9-Shear stress sampled on intennittency through transition (Case 2). 
A * airnnar,. # - turbulent,+ * transition. Station 3a refers to data taken at 
station 3, but at a reduced free-stream velotity. 
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Fig. 3.2.10-Turbulent heat flux sampled on intermittency through transition 
(Case 2). A- laminar, • • turbulent, + - transition. The intermittency at 
station 2 is 9.5%. The other profiles are fully turbulent. 
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Fig. 3.2.11— Turbulent Prandtl number sampled on intermittency through 
transition (Case 2). A* laminar, # - turbulent, +* transition. The 
intermittency at station 2 is 9.5%. The other profiles are fully turbulent. 
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3.3). Case 3 — Flat-wall, TI=8.3% 


The measurements described in this section were taken with the 
turbulence generating jet grid (described in Chapter 2) in place. The free- 
stream turbulence generated, though, was so high that it was difficult to 
distinguish between laminar and turbulent flow in the transition region from 
the hot-wire traces. The situation is illustrated on Fig. 3.3.1 where a hot-wire 
voltage trace in the transition region and it's first derivative are presented. 
The derivative of the signal seems to fall into two distinct regimes of 
fluctuation amplitude, suggesting that the flow is transitional, but this is 
difficult to verify from the direct hot-wire trace. If the velocity is increased 
(thereby increasing the local Reynolds number), the large-amplitude- 
fluctuation regime expands at the expense of the low-amplitude-fluctuation 
regime, further suggesting a transitional flow. The two regimes are hot very 
distinct, however, making separation of the flow into laminar and turbulent 
regimes difficult (a massive number of dropouts occurs). For this reason, no 
intermittency based processing was performed for this case. 

Free-stream Turbulence Intensity and Spectra . The power spectral 
density measured in the free-stream at St. 1 is shown on Fig. 3.3.2. The 
spectrum is seen to be clean, with no significant spikes, in contrast to the two 
lower TI cases. The high free-stream turbulence generated by the jet grid 
apparently overwhelms any tunnel unsteadiness and/or electronic noise. 

The free-stream turbulence intensity at the tunnel centerline was 
found to be quite isotropic, as seen from Fig. 3.3.3, and decays from 8.3% at St. 
1 to 5.9% at St. 4. There was a significant variation in v' in the cross-stream 
direction, however, with v' changing by as much as 45% from just outside the 




boundary layer to the tunnel centerline. This is thought to be due to the 
particular jet-grid geometry used. 

Stanton Number . The wall Stanton number variation is shown on 
Fig. 3.3.4. The data is seen to deviate from the turbulent correlation for 
Re x <lxl05, indicating a short "transition" region, consistent with the 
intermittent hot-wire signal discussed above. The laminar region (if one 
exists) is too short to be measured. 

In contrast to the data of Blair (1983) and Simonich and Bradshaw 
(1978), the increase in free-stream turbulence is not seen to augment the heat 
transfer. This may be due to the low Reynolds number, however, as 
suggested by Simonich and. Bradshaw (1978). An energy balance is shown in 
Fig. 3.3.5. The heat transfer closure is excellent. 

Mean Velocity Profiles . Profiles of mean velocity plotted in wall 
coordinates are shown on Fig. 3.3.6. The profiles corresponding to stations 2, 
3, and 4 are seen to agree very well with the log-linear law. As shown, the 
log-linear region expands with Re x . A notable feature of the profiles is the 
absence of a wake. This is due to the high free-stream turbulence level. The 
profile corresponding to station 1 is seen to be transitional. As there was 
some difficulty in assigning an appropriate Cf value to this profile due to the 
absence of a log-linear region, the profile shown may not be precisely correct. 
Determining the local skin friction for this profile using a momentum 
balance was not possible since no data upstream of St. 1 was taken. 

Plots of the skin friction and shape factor (H) are shown on Figs. 3.3.7 
and 3.3.8, respectively. The skin friction is seen to drop monotonically with 
Re x . The shape factor was measured to be 1.71 at St. 1, indicating a 
transitional boundary layer at this location, but one which has nearly reached 
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a fully-turbulent shape. The shape factor, H, equal to 1.4, the accepted value 
for turbulent boundary layers at stations 2, 3, and 4. 

Velocity Fluctuation . Fluctuations of the streamwise turbulence 
intensity are shown on Fig. 3.3.9. The peak in the profile at St. 1 is relatively 
broad, and is indicative of a laminar or transitional flow. The peaks for the 
profiles at stations 2, 3, and 4 are much sharper, with the peaks dropping 
monotonically with Re x and in step with the drop in free-stream turbulence 
level with successive downstream positions. 

Mean Temperature Profile . Mean temperature profiles measured 
using a thermocouple probe and normalized on wall coordinates are shown 
on Fig. 3.3.10. The profiles at stations 2, 3, and 4 show log-linear regions. As 
was seen for the corresponding velocity profiles, the extent of the log-linear 
region increases with Re x . No wake is seen. Turbulent Prandtl numbers 
deduced from the profiles were consistently near unity. The temperature 
profile at St. 1 does not seem to possess a log-linear region, consistent with the 
transitional nature of the flow at this station. This profile was obtained using 
the uncertain value of Cf obtained from the mean velocity profile. 

Shear Stress Profiles . Profiles if u'v’ are shown on Fig. 3.3.11. The 
profile at St. 1 with it's broad peak is indicative of a transitional flow profile, 
while the profiles at station 2, 3, and 4 have a turbulent flow shape. The near- 
wall peak is seen, to decrease with Re x . 

Turbulent Prandtl Number Measurements . Profiles of the turbulent 
Prandtl number (Prt) measured using a special triple-wire probe are presented 
on Fig. 3.3.12. Prt values are greatly increased above unity for the early 
turbulent boundary layer (Stations 2 and 3), but are seen to decay to nominally 
unity by station 4. This suggests that the momentum boundary layer 
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establishes itself more quickly than does the thermal boundary layer in the 
early turbulent flow, resulting in higher values of the eddy diffusivity of 
momentum relative to the eddy diffusivity of heat. This view is supported by 
the mean velocity and temperature profiles, where similar velocity profiles 
were seen at stations 3 and 4 while the temperature profiles were still 
evolving. There is no reason to expect the momentum and thermal 
boundary layers to develop at the same rate, since the boundary conditions 
are different. The momentum boundary layer sees a non-zero fluctuation (a 
non-zero u') in the free-stream whereas t' in the free-stream must equal zero. 
It is postulated that having eddies present in the free-stream enables the 
momentum boundary layer to respond more quickly than if it had to grow by 
turbulent diffusion alone. Creating a temperature fluctuation in the free- 
stream (possibly by injecting heated air through the grid) may cause the 
thermal ard momentum boundary layers to grow at comparable rates. Prt 
values are in the vicinity of unity by station 4, indicating that the momentum 
and thermal boundary layers have grown to comparable thicknesses. 


85 



• 1 i 1 1 ' 1 ' 

2 


> VL 
1—1 

in 


o 

> 


i 




-«-■ _ I i I 1 L 


_L 


Hot-wire signal 


1 st derivative 


.01638 .03276 .04514 .06552 .0819 

TIME [S3 


a). Near free-streaxn 



Hot-wire signal 


— < — - I . i _ I — Li i i L 


1 st derivative 


.0l63e .03276 .04914 .06552 .0819 

TIME CS3 


b). Near-wall flow 


Fig. 3.3.1-Hot-wire trace and its derivative in the boundary layer (Case 3). 
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Fig. 3.3.2-Power spectral density distribution of streamwise velocity (Case 3). 
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Fig 3.3.4— Stanton number variation along test-wall (Case 3). 
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Fig. 3.3.5-Energy balance along wall (Case 3). 
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Fig. 3.3.6— Mean velocity profiles normalized on wall coordinates (Case 3). 
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Fig. 3.3.9-Turbulence intensity profiles along the test-wall (Case 3). 
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Fig. 3.3.12-TurbuIent Prandtl number profiles along the test- wall (Case 3). 
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3.4). Case 4 — Concave-wall. TI=0.60% 

The data for this case were taken with the flexible heated wall bent to a 
radius of curvature of 97 cm. No turbulence generating grid was used. The 
outer flexible wall was adjusted such that the static pressure at the test wall 
was uniform within ±3% of the upstream static pressure. A comparison of 
the measured free-stream velocity variation within the curve at stations 3 
and 4 with the theoretical velocity distribution as computed by potential 
theory is shown on Fig. 3.4.1. The theoretical velocity distribution was 
computed by fitting the equation 


Up. 

u ^ )= — T 

R (2.6) 

through two of 'he measured points. It is seen that although there is a slight 
discrepancy between the measured and theoretical profiles (the source of the 
discrepancy will be discussed in section 3.5— the curved wall, high turbulence 
case— where the discrepancy is much more severe), the two agree quite well. 
Because of the discrepancy between the two, however, the calculation of 
boundary layer thicknesses is not precise. For this reason, all distances 
normal to the wall in this section will be normalized on the wall curvature 
(R). This is done also to be consistent with the high-disturbance cases where 
integral thicknesses and the 99.5% boundary layer thickness cannot be 
evaluated (as will later be discussed). 

A plot of the hot-wire signal in the transitional flow along with the 
corresponding intermittency signal is shown on Fig. 3.4.2. A large number of 
dropouts is observed, indicating that it is difficult for the circuit (and also the 
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observer) to determine which portion of the hot-wire signal represents 
turbulent-like flow and which portion represents laminar-like flow. Close 
inspection of the hot-wire signal revealed that it contained a relatively small 
amplitude, high frequency signal superimposed on a large amplitude, low 
frequency signal. It was difficult to decide whether to classify this as a 
turbulent flow or an unstable laminar flow. Furthermore, the derivatives in 
the signal were of comparable magnitudes, making it difficult for the circuit 
to decide whether the flow was laminar or turbulent. For these reasons, 
conditional sampling was not performed for this case. Since transition occurs 
so rapidly, it was possible to adjust the free-stream velocity such that 
transition occurred between stations (stations 2' and 3). This enabled data to be 
taken in the late laminar (station 2) and early turbulent (station 3) flows, 
completely bypassing the transition region. 

A photograph of the liquid crystal on the test wall (with wall heating— 
Up W ~ 17.2 m/s) is shown on Fig. 3.4.3. Transition occurs (for this wall heat 
flux) close to station 2 via, it appears, a breakdown of the vortex structure. 
The spanwise variation implies that the vortices break down independently 
of each other. This is consistent with the observations of Swearingen (1985), 
who used smoke to visualize the vortex breakdown. The transition pattern 
seen here is quite unlike that observed on the flat-wall, in which little 
spanwise irregularity in the transition pattern was observed. The spanwise 
temperature variation in the laminar flow is relatively stable in time and 
stationary in space, implying that the vortices are also stable. If the vortices 
do move, and it is quite possible they do, they move at a frequency larger than 
the frequency response of the liquid crystal (~1 Hz) with an amplitude smaller 
than the vortex wavelength. Spanwise variations in temperature are also 
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seen further downstream (stations 3 to 6), implying that a relatively stable 
vortex structure exists in the turbulent flow as well. This is at odds with the 
results of Simonich and Moffatt (1982) who concluded that a stable three- 
dimensional vortex structure did not exist bn the concave wall when a 
turbulent boundary layer was introduced to concave curvature. Barlow and 
Johnston (1988), in a similar study, found that vortices existed, but were not 
fixed in time and space, the vortices appearing and disappearing randomly 
across the span. Only when vortex generators were used upstream of the 
curve was the vortex position fixed. It appears that in the present study, the 
stable vortex structure in the laminar flow serves the equivalent function to 
that of the vortex generators in Barlow and Johnston's (1988) study, in that 
they provide preferred spanwise locations for the vortices. The advantage of 
the present study is that this "locking" is done "naturally". Lastly, the vortex 
wavelength is seen to be quite irregular across the span in both the laminar 
and turbulent flows. It was thought that the geometry of the step at the 
leading edge of the liquid crystal sheet might have an effect on the vortex 
spacing. Additional effort was expended on smoothing the joint by sanding 
and a picture of the vortex pattern was re-taken. Comparison of the vortex 
pattern before and after this sanding revealed no difference between the 
patterns. It was therefore concluded that the joint had little influence on the 
vortex spacing. The parameter controlling the vortex spacing is probably the 
last screen upstream of the nozzle, as was concluded by other researchers (e.g., 
Swearingen-1985, and Bippes— 1978). 

Wall heating was found to affect the transition location , as was seen in 
Case 1. A plot of the intermittency vs. time for wall heating and wall cooling 
is shown on Fig. 3.4.4. The intermittency is seen to vary from 40% with no 
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wall heating to approximately 90% with wall heating. As the wall cools, the 
intermittency returns to 40%. As in Case 1, this precludes the use of velocity 
profiles measured in an unheated flow to reduce the temperature profile data 
(such as the computation of enthalpy thickness), at least in the laminar and 
transitional flow cases. It is felt that this problem does not exist in the post- 
transitional profiles, and so the data at these stations (3, 4 and 5) will therefore 
be presented. 

Wall heating destabilizes the flow in two ways. First, wall heating 
increases the viscosity near the wall, leading to inflectional velocity profiles as 
was discussed in Case 1. Secondly, the fluid density close to the wall 
decreases, causing the heated fluid to move away from the wall due to the 
centrifugal forces in the channel. An estimate of the effect of changes in fluid 
density on flow stability can be calculated using the results of Lin, Kamotani 
and Ostrach (1982), who investigated the effects of buoyancy on Gortler 
vortices in a concave-curved channel heated from below. They found that for 


where Gr is the Grashof number based on gravitational acceleration, 
buoyancy forces dominated the centrifugal forces and an appreciable increase 
in the vortex amplitude occurred. The above parameter, the ratio of 
buoyancy to centrifugal forces, computed for the present case using the 
centrifugal acceleration in place of the gravitational acceleration (the 
centrifugal acceleration was over 30 times greater than the gravitational 
acceleration) was found to be 


101 



implying that the buoyancy term has little effect on the observed instability. 
The destabilizing effect of heating on the flow thus seems to be caused mainly 
by the increase in fluid viscosity near the wall. 

Another potential problem with the temperature data occurs during 
the normalization to wall coordinates, T + vs. y + , where the local wall heat 
flux and temperature are required. Unfortunately, the wall geometry 
prevents the accurate determination of these quantities when the flow is 
strongly three-dimensional. As shown on Fig. 2.4, the thermocouples are 
embedded behind a ~1 mm (40 mil) lexan/liquid crystal composite. This 
composite tends to smear out temporal and spatial variations in wall 
temperature, resulting in much smaller variations at the thermocouple 
locations. The temperatures recorded by the thermocouples are thus averages 
of the local surface temperature variation, and do not represent the true 
surface temperature. A similar problem occurs with the heat flux due to the 
conduction of heat within the composite. The heat flux is also a measure of 
the average and not the local' value. The problem discussed above is 
significant only when the order of the non-uniformity (the vortex spacing in 
the present case) is comparable to or smaller than the composite thickness. 
The vortex wavelength in the late laminar flow was found to be ~4 mm, 
giving a disturbance wavelength to composite thickness of ~4. Though this is 
encouragingly large, some smearing of the temperature and wall heat flux 
variation is expected. Since the vortex wavelength increases from 10-25 mm 
after transition, however, the problem becomes much less severe and 
accurate values of the local wall temperature and heat flux are obtainable. 



Due to the problems encountered in the laminar flow, only post-transitional 
profiles will be presented. 

The embedded thermocouples at the centerline of the test wall were 
found to lie almost directly beneath a downwash. The wall temperatures 
beneath the downwash could thus be found. To find the wall temperatures at 
the upwash, two stick-on foil thermocouples 0.0127 mm (0.5 mils) thick 
(made by Rdf Corporation), configured to give the temperature difference 
between the two junctions, were taped onto the wall at the upwash and 
downwash locations. Knowing the wall temperature difference between the 
upwash and downwash and the wall temperature at the downwash, the wall 
temperature at the upwash could be calculated. 

Mean and Fluctuating Velocity . Results Of various spanwise traverses 
of the hot-wire at constant y-distances from the wall for stations 2 to 5 are 
shown on Fig. 3.4.5. The spanwise variation of the mean (Fig. 3.4.5) and 
fluctuating velocities (not shown) in the laminar flow (St. 2) is especially 
pronounced, with a peak in velocity corresponding to a dip in fluctuating 
velocity, and vice versa, illustrating the unstable nature of the inflection 
point velocity profile in the upwash. Oscilloscope traces of the hot-wire 
output at a constant y-distance from the wall in the upwash and downwash, 
shown on Fig. 3.4.6, illustrates the difference in fluctuation level between the 
two flows. It is further seen from Fig. 3.4.5 that the vortex spacing in the late 
laminar flow (station 2) is quite irregular, as was also seen in the liquid crystal 
visualization, and that the upwash and downwash do not line up from one y- 
distance to another, indicating a tilted vortex structure. Also, close inspection 
of the fluctuating velocity revealed a double peak in u' within the vortex, 
consistent with the observations of Swearingen (1985). 
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The location of the state of the flow at station 2 on the Gortler map, 
shown on Fig. 3.4.7, was determined from calculated values of G and k82- 
The momentum thickness corresponding to the average between the upwash 
and downwash values was used in the calculations. The location on the map 
is well within the unstable range. Also plotted on Fig. 3.4.7 are lines of non- 
dimensional wavelength. A, given by 

A _ ttoo R m 3/2 

A “ v 1 r ) (3.2) 

which were calculated from stability theory. From the spanwise 
measurements, an estimate of the wavelength (X) was made from which the 
nondimensional wavelength. A, was calculated to be 281. This is in good 
agreement with the map values. 

Shown on Fig. 3.4.8 is a plot of the transition Reynolds number vs. a 
non-dimensional radius (R/82) from Liepmann (1943). The flat plate 
transition Reynolds number is shown by a horizontal line, while Liepmann's 
criterion is shown by the sloped line marked N=240 (N is the Gortler number, 
G, calculated using the boundary layer thickness in place of the momentum 
thickness). Note that for smaller radii, i.e., smaller R/82,tr/ the data fit 
Liepmann's stability criterion. However, for larger radii, the data approach 
the flat plate transition criteria. It can therefore be concluded that transition 
to turbulence on highly concave surfaces is dominated by the Gortler vortices. 
Also marked on Fig. 3.4.8 is the state of the present flow. It is seen to fall near 
Liepmann's stability criterion, implying that the Gortler vortices dominate 
the transition process. 
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The difference in the span wise profiles between the upwash and 
downwash (Fig. 3.4.5) becomes much less pronounced after transition, due to 
the increased mixing in the boundary layer. The vortex wavelength becomes 
much larger, and no double peak in u', as was observed in the laminar flow, 
is seen. The crests and troughs also align. 

The growth of the boundary layer and the effect of the vortices on the 
mean velocity profiles are illustrated on Fig. 3.4.9. Large differences between 
the upwash and downwash profiles are seen in the late-laminar flow (station 
2) with a clearly inflectional upwash velocity profile. The profiles become 
quite similar just after transition, then again deviate at the later stations. An 
explanation for this remarkable behavior wiH, be discussed later. The mean 
velocity normalized on wall coordinates at the upwash and downwash 
locations near the tunnel centerline are shown on Fig. 3.4.10. The upwash 
profile at station 2 is distinctly inflectional and a dramatic difference between 
the upwash and downwash profiles is seen. Again, much of this difference 
disappears after transition. Although it appears that the upwash profiles 
have a fuller shape than the downwash profiles, this is an artifact of the 
normalization (the skin frictions for the upwash are smaller than for the 
downwash, resulting in higher u + values in the upwash). 

Plots of the skin friction and shape factor vs. Re x are shown on Figures 
3.4.11 and 3.4.12, respectively. In both plots, large differences in values in the 
late-laminar region, followed by closer values in the turbulent flow are seen, 
supporting the trends seen in both the spanwise and normal profiles of 
streamwise velocity. It is interesting to observe, however, that the upwash 
and downwash values approach one another just downstream of transition, 
then increasingly deviate downstream, implying that the vortices again 
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become stronger. This led the authors to speculate that this behavior may be 
reflected in the "turbulent" Gortler number (Gt) (the Gortler number with the 
eddy viscosity replacing the molecular viscosity). This value decreases to a 
low value immediately after transition due to the sudden increase in eddy 
viscosity, then slowly increases in value as the turbulent boundary layer 
grows. Measurements of G t using a cross-wire probe to measure the shear 
stress and the mean velocity gradient to obtain the eddy viscosity (v e ) are 


shown below. 


Station Sox'! 000 fml Ve [m^ /s] v\m 2 M Gj 


2 (upwash) .530 

2 (downwash) .164 

3 (upwash) .998 

3 (downwash) 1.124 

4 (upwash) 1.167 

4 (downwash) 1.820 

5 (upwash) 1.898 

5 (downwash) 2.718 


.000016 13.06 

.000016 2.230 

0.000425 1-288. 

0.00043 1-522 

0.000894 0.592 

0.000882 1-531 

0.001046 1-545 

0.001329 1-814 


The eddy viscosity used is the average value across the boundary layer. It was 
computed by numerically integrating the values then dividing by the 
boundary layer thickness. The measurements show that Gt does decrease to 
~1.3 immediately after transition, but remains fairly constant thereafter. G t 
values computed using the maximum value of eddy viscosity within the 
boundary layer decrease the magnitude of Gt, but show a similar trend. 

Profiles of u’ are shown on Fig. 3.4.13. The large variation in u in the 
late laminar flow between the upwash and downwash is seen to completely 
disappear after transition (station 3), where similar profiles are seen. The 
differences re-emerge at the later stations (stations 4 and 5) in the outer part of 
the flow, consistent with the observed trends in mean velocity, Cf and H. The 
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near-wall peak in u' is consistently near 11% for all the post-transitional 
profiles indicating a near-wall curved-asymptotic situation. 

Shear Stress Profiles . Profiles of the shear stress u'v' are shown on Fig. 
3.4.14. The data at station 2 is not considered reliable since the spacing 
between the wires of the X-wire probe is comparable to the vortex 
wavelength. When the probe was placed in the boundary layer and traversed 
in the spanwise direction, it was clearly seen that both wires were not in a 
upwash or downwash simultaneously. The data at this station was taken by 
centering the probe at the upwash or downwash as best the author could. 

The shear stresses at the wall plotted on these figures are not measured 
values of u'v' but computed values obtained from the skin friction values 
deduced from the mean velocity profiles. Generally good agreement between 
the wall values and the profile values are seen. A reversal in shear stress in 
the upwash at station 2 is observed due to the inflection point velocity profile. 
The shear stresses in the upwash are greater than the shear stresses in the 
downwash in the turbulent flow, at odds with what the skin friction values 
would suggest (Cf in the upwash is smaller than in the downwash). A 
distinction must be made between the riear-wall flow and the wake flows, 
however. The wake flow is significantly affected by the large vortical motion 
which convects turbulent fluid towards the upwash. The near-wall flow is 
less affected, depending more on the local velocity profile. 

Stanton Number . The Stanton number variation along the wall under 
an upwash is plotted on Fig. 3.4.15 along with the corresponding flat-plate 
results. Concave curvature is seen to be highly destabilizing, causing 
transition to occur about five times earlier than on the flat plate. It was noted 
that the transition start, path and length varied depending on whether the 
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centerline thermocouples were beneath a downwash or an upwash, however. 
An example is given on Fig. 3.4.16 where Stanton number measurements at 
two free-stream velocities were made. The change in the boundary layer 
thickness causes the vortex wavelength to change, enabling the embedded 
thermocouples along the centerline to lie under an upwash for the 6.74 m/s 
case or a downwash for the 17.2 m/s case. Transition is seen to occur over a 
shorter length under an upwash than under a downwash. Stanton number 
values for the upwash locations in the laminar flow are seen to lie below the 
laminar correlation. It is not known whether the downwash values lie 
slightly above the correlation due to unheated starting length effects. 

Mean and Fluctuating Temperature Profiles . Mean temperature 
profiles normalized on wall coordinates are plotted on Fig. 3.4.17. The 
striking feature of these profiles is their deviation from the thermal law-of- 
the-wall. Pauley and Eaton (1988) found that for a pair of spanwise-separated 
vortices embedded in a boundary layer where the common flow of the 
vortices is toward the wall, a strong increase in T + was observed for all 
locations within the vortex except in the upwash at the edge of the vortex pair 
where a decrease in T + was observed. The increase in T + values in both the 
present case and in the study by Pauley and Eaton (1988) is thought to be due 
to dilution of the heated boundary layer flow by the free-stream flow, 
resulting in an overall lowering of the temperature in the boundary layer. 
The difference between the wall temperature and the temperatures in the 
boundary layer increases, resulting in higher values of T + . The present 
results and those of Pauley and Eaton (1988) indicate that the thermal law-of- 
the-wall is not valid in boundary layers which have such large, streamwise 
vortices embedded within them. 
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Profiles of the fluctuating temperature are shown for the post- 
transitional stations on Fig. 3.4.18. The data was obtained using the triple- 
wire probe described earlier. Values of t’ in the upwash and downwash are 
similar for station 3, then deviate for stations 4 and 5, similar to the behavior 
observed for the corresponding u' profiles. Peak values of t' are —12% of the 
wall to free-stream temperature difference for all stations. The temperature 
fluctuations in the upwash are greater than those in the downwash, reflecting 
the relatively unstable nature of the flow and the convection of heated flow 
(and thus larger fluctuation magnitudes) toward the upwash. 

Turbulent Heat Flux and Turbulent Prandtl Number . Profiles of the 
streamwise and normal cross-stream heat fluxes normalized on the wall heat 
flux are shown on Figs. 3.4.19 and 3.4.20, respectively. The streamwise heat 
flux is roughly twice the cross-stream heat flux for all profiles at both the 
upwash and downwash locations. The cross-stream heat flux profiles 
approach unity near the wall, as expected. The cross-stream diffusion of heat 
is greater in the upwash than in the downwash, similar to the behavior 
observed in the shear stress profiles. Again, this occurs due to span wise 
convection of heat to the upwash locations. Values of v't' in the upwash are 
greater than those in the downwash even though Stanton numbers in the 
upwash are lower than those in the downwash. 

Profiles of the turbulent Prandtl number deduced from u'v', v't', and 
the gradients in velocity and temperature are shown on Fig. 3.4.21. All 
profiles are seen to be near unity in the vicinity of the wall, indicating no 
gross violation of Reynolds analogy. The data in the outer part of the 
boundary layer is not reliable due to the weak values of u'v' and v't' and the 
very shallow gradients in velocity and temperature. 
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Fig. 3.4.2 Hot-wire signal and output of intermittency drcuit (Case 4). 
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Fig. 3.4.3— Picture of heated test-wall showing spanwise variation in wall 
temperature (Case 4). 
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Fig. 3.4.4— Variation of intermittency vs. time for wall heating and wall 
cooling (Case 4). 
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Fig. 3.4.5-Spanwise variation in velocity at various normal distance from the 
wall ffacA dV 
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Fig. 3.4.6— Oscilloscope traces at upwash and downwash, station 2 (Case 4). 
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Fig. 3.4.7— Location on Gortler map of the late laminar (station 2) flow (Case 4). 



Fig. 3.4.8— Location on Liepmann's (1943) transition diagram of the late 
laminar (station 2) flow (Case 4). 


116 





117 






Fig 34.10-Velodty profiles at the upwish and downwash normalized on 
wall coordinates (Case 4). 
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Fig. 3.4.11-Skin friction variation along the wall at upwash and downwash 
locations (Case 4). 
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Fig. 3.4.13-Turbulence intensity profiles along the wall at upwash and 
downwash locations (Case 4). 
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Fig. 3.4.14-Shear stress profiles along the wall at upwash and downwash 
locations (Case 4). 
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Fig. 3.4.15— Stanton number variation vs. Rex— comparison with the low TI 
flat-wall case. 
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Fig. 3.4.16-Stanton number variation vs. Rex for two free-stream velocities 
(Case 4). 
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Fig. 3.4.18-Fluctuating temperature profiles along the wall at upwash and 
downwash locations (Case 4). 


126 




3 



0 0.02 0.04 

y/R 


Fig. 3.4.19-Profiles of the stream wise heat flux along the wall at upwash and 
downwash locations (Case 4). 
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Fig. 3.4.20-Profiles of the cross-stream heat flux along the wall at upwash and 
downwash locations (Case 4). 
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Fig. 3.4.21-Profiles of the turbulent Prandtl number at upwash and 
downwash locations (Case 4). 
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3.5). Case 5 — Concave-wall. TI=8.3% 

This case deals with the effects of concave curvature on transitional 
boundary layers under high free-stream turbulence conditions. The free- 
stream turbulence at the test section entrance, generated using the bi-plane 
grid generator described earlier, was measured to be 8.6%, similar to the 
corresponding flat- wall, high turbulence intensity case (Case 3). The PSD 
distribution was smooth, with no significant spikes over the range from 10 
Hz to 10 kHz. Velocity and turbulence intensity profiles at the end of the 
turbulence establishment chamber and just before the test section entrance 
(see Fig. 3.5.1) showed a mean velocity variation across the span of 2% and a 
variation in turbulence intensity of 6%. The autocorrelation, given by 


p(t) = 


u'(t)u’(t + t) 
u'2(t) 


can be used to find two turbulence scales. The area under the autocorrelation 
curve, called the integral scale, represents the average size of the turbulent 
eddies. The Taylor microscale, which is related to the turbulent dissipation, is 
determined from the curvature of the autocorrelation curve at the origin 
(t= 0). A description of this measurement, written by Mr. Steve Russ, is given 
as an appendix. The integral scale and Taylor microscale at the test section 
entrance were measured to be 3.3 cm and 0.61 cm respectively at the center of 
the channel. 

Perhaps the most startling find of this case was the phenomenon of 
cross-stream transport of momentum within a flow that was thought to be a 
potential core. Apparently, the combination of a high free-stream turbulence 
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intensity superimposed on a free-stream velocity gradient (due to 
conservation of angular momentum within the curve) causes a transport of 
momentum within the "potential core", as manifested in high levels of shear 
stress (u'v') . As a result, there exists a velocity profile in the core that is 
flatter than the profile predicted by potential theory (given by 
U(y)r(y)=constant) due to the increased mixing. Mean velocity profiles, 
measured across the test section normal to the test wall at each station (Fig. 
3.5.2) increasingly deviate from the potential flow distribution with 
downstream distance. The deviation is seen as early as station 2, where one 
would expect the boundary layer to still be thin. The profiles are seen to be 
flatter than the potential flow profile. It is hypothesized that this is due to a 
large momentum transport in the "core" flow. This seems to be supported by 
the shear stress profiles (Fig. 3.5.3) where large values of shear stress are seen 
even at the channel centerline. It seems that the high turbulence intensity in 
the core, when superimposed on a velocity gradient, causes transport of 
momentum from the flow near the convex wall (higher velocities) towards 
the flow near the concave wall (lower velocities). The production term in the 
shear stress budget equation, given by 


P=v^-k“' 2 -v' 2 )# 

shows that the production of shear stress can be positive for non-zero free- 
stream turbulence when streamline curvature is present (note that for 
concave curvature, R<0). The parameters thought to govern this 
phenomenon are the free-stream mean velocity gradient and turbulence 
level, the wall curvature and the channel width. 
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A consequence of this was that the usual normalizing techniques were 
not applicable since neither a potential velocity at the wall nor a boundary 
layer thickness (there was no boundary layer) existed. The quantity selected to 
normalize velocities was the potential velocity at the wall as determined by 
the upstream total pressure and the radius of curvature of the wall, i.e., a 
velocity profile in the free-stream was inferred using the measured total 
pressure upstream of the curve, and the velocity at the wall was the quantity 
chosen. Normal distances from the wall were normalized on the wall 
curvature, R. Momentum balances were not possible in this case, due to the 
cross-transport of momentum, unless detailed profiles were taken near both 
the convex and concave walls for evaluation of wall shear at both locations. 
The present facility did not allow such measurement. Energy balances; 
however, could be made if velocity and temperature profiles were measured 
beyond the thermal 'boundary layer" thickness. 

Another interesting result is that no evidence of streamwise vortices 
were seen. The turbulent Gortler number (Gt ) could not be calculated for this 
case as no momentum thickness could be obtained. The values of the eddy 
viscosity in the high turbulence intensity case flow were, however, anywhere 
from 10 to 100 times that of the lower turbulence intensity case, suggesting 
that if Gt could have been calculated, they would have been much lower than 
the values computed for the low turbulence intensity case. Possibly, Gt would 
have been reduced to the stable region on the Gortler map, and no vortices 
would would be expected. Alternatively, it could be that vortices exist, but 
that they meander or appear and disappear randomly in time and space at a 
rate faster than the frequency response of the liquid crystal. Clearly, the flow 
appears two-dimensional. 
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Mean and Fluctuating Velocity . Measurements of the mean velocity 
normalized on wall coordinates are shown on Fig. 3.5.4. As in the 
corresponding flat-wall case, no wake, and a very short log-linear region is 
observed. Both concave curvature and high free-stream turbulence reduce 
the size of the wake. Both effects are visible in these profiles. Station 1 is 
thought to be transitional. Skin friction (Cf) values deduced from the law-of- 
the-wall showed a monotonic decrease with Re x . 

The variation in turbulence intensity profiles along the test wall is 
shown on Fig. 3.5.5. The near-wall peak decays slowly with distance. The 
turbulence intensity in the outer part of the flow decays rapidly at first, and 
essentially stops decaying beyond station 3. Baines and Peterson (1951), who 
studied the decay of turbulence behind screens, found that the decay of free- 
stream turbulence behind lattice type grids was governed by the equation 



where x is the distance from the grid generator and b is the bar width. The 
above equation was used to find the effective bar width for the present flow 
(note that cylinders instead of bars were used in the turbulence generator) by 
solving for the bar width which yielded the measured u7u„ at station 1. The 
effective bar width was found to be 2.78 cm compared to the actual cylinder 
diameter of 4.45 cm. The equation was then used to predict what the 
turbulence decay rate would have been if the channel had been straight. 
Results are plotted on Fig. 3.5.6. It is seen that the turbulence in the straight 
section would have continued to decay if it were not for the cross-stream 
transport of momentum in the curved-wall flow. In contrast to the curved- 
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wall flow, the high free-stream turbulence levels in the flat-wall flow (Fig. 
3.3.9) showed a continuous decay. 

Stanton Numbers . Stanton numbers along the wall centerline are 
shown on Fig. 3.5.7. A line representing the data obtained for the high 
turbulence intensity, flat-wall case is also shown. Concave curvature is seen 
to increase the heat transfer from the wall significantly, suggesting increased 
mixing of the near-wall flow. Vortex motion may be responsible for this 
increase, although the high eddy viscosity discourages the formation of 
coherent vortices. If vortex motion does exist, it must be quite disorganized 
as no evidence of a spanwise variation in heat transfer was observed on the 
liquid crystal. Alternatively, the concave curvature may lead to increased 
instability and more frequent turbulence bursting, but not coherent cellular 
structures. Thus, this concave- wall flow is two-dimensional. An energy 
balance for this flow (Fig. 3.5.8) shows exce’lent closure. 

Mean and Fluctuating Temperature Profiles . Mean temperature 
profiles normalized on wall coordinates are shown on Fig. 3.5.9. As in the 
lower turbulence intensity case, a discrepancy between the measured values 
of T + and the thermal law of the wall is seen. The discrepancy is much 
smaller in the present case, however, illustrating the effects of increased 
turbulence intensity. 

The variation in fluctuating temperature, measured using the triple- 
wire, is shown on Fig. 3.5.10. Unlike the turbulence intensity profiles, the t' 
profiles are seen to approach zero in the outer portion of the flow since the 
"core" flow is isothermal. This illustrates a fundamental difference between 
the momentum and heat transfer processes. The boundary conditions for the 
two processes are different (similarity in boundary conditions in the high 
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turbulence intensity case could be achieved if the outer wall were heated and 
heated air were injected into the free-stream). Due to this non-similarity, the 
turbulent Prandtl number is not expected to equal unity. 

Turbulent Heat Flux and Turbulent Prandtl Numbers . Profiles of the 
stream wise and cross-stream transport of heat are shown on Figs. 3.5.11 and 
3.5.12, respectively. Both profiles show an evolution with downstream 
distance as heat diffuses away from the wall. The v't' profiles approach unity 
near the wall, as expected. In contrast to the shear stress profiles, which 
remained high across the test section, the turbulent heat flux profiles 
approach zero in the outer part of the flow. The difference in boundary 
conditions between the heat and momentum transfer processes is again 
illustrated. 

Turbulent Prandtl numbers deduced from the triple-wire 
measurements are shown on Fig. c.5.13. All the near- wall values are seen to 
be slightly higher than unity. This increase in Pr t is not surprising given the 
difference in boundary conditions discussed above. 
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Fig. 3.5.1-Profiles of the mean and fluctuating velocities across the test section 
at the end of the turbulence establishment chamber (Case 5). 
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Fig. 3.5.5-Turbulence intensity profiles along the test wall (Case 5 ). 
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Fig. 3.5.6— Comparison of the turbulence decay rate in the curved test section 
with the predicted decay rate of Baines and Peterson (1951) in straight 
channels. 
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Fig. 3.5.7-Comparison of the Stanton number variation along test wall 
between curved and straight wall cases. 


142 


1E+04 1E+05 IE+06 1E-H37 

Rex 





[unuj 

Fig. 3.5.8— Energy balance for Case 5. 
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Rg. 3.5.10— Profiles of the fluctuating temperature along the test wall (Case 5). 
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Fig. 3.5.11-Profiles of the stream wise heat flux along the test wall (Case 5). 
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Fig. 3.5.12-Profiles of the cross-stream heat flux along the test wall (Case 5). 
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Fig. 3.5.13-Profiles of the turbulent Prandtl number along the test wall (Case 

5 ). 
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CHAPTERS 


Conclusions 

The effects of free-stream turbulence and concave curvature on transitional 
boundary layers were studied. The main conclusions of this study are : 

1) . The flat-plate transitional boundary layer cannot be thought of as being a 
simple composite of a Blasius and a fully-turbulent flow. Transition 
modelling based on the intermittency function weighting of pure laminar 
and turbulent flows may be in error. 

2) . Conditional sampling of turbulence quantities on the intermittency 
function must be made during transition. Measurements of time-averaged 
quantities may not give an accurate view of the transition process. 

3) . The turbulent Prandtl numbers in the turbulent core region of the 
transitional flow are somewhat smaller than unity. 

4) . The existence of stable vortices on the concave-curved wall in both the 
laminar and turbulent boundary layers was established for low free-stream 
turbulence intensities. No coherent vortices were found for the higher 
turbulence intensity case. 

5) . Concave curvature destabilizes the flow, causing transition to occur earlier 
than on the flat-wall. This is a confirmation of earlier findings. 

6) . No gross violation of Reynolds analogy was found for the post- 
transitional profiles in both the low and high turbulence intensity cases 
although small deviations from an exact analogy were noted. 

7) . High levels of free-stream turbulence superimposed on a free-stream 
velocity gradient were found to cause a cross-stream transport of momentum 
within the "potential core" of the flow. 
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Al). Measuremen t of Emissivitv 

The emissive power of a surface is defined as the radiant energy 
emitted by a surface per unit time and unit surface area. A blackbody absorbs 
all the radiation incident upon it and is the most efficient emitter of 
radiation. The emissive power of a blackbody is given by the Stephan- 

Boltzmann law 


where 


qb-tfl 4 

q b = emissive power [W/m 2 ] 
a = Stephan-Boltzmann constant 
T = absolute temperature [K] 


(Al.l) 


Most surfaces emit less radiation than that emitted by a blackbody at the same 
temperature. The ratio of the emissive power of a non-blackbody to that of a 
blackbody is called the emissivity of the surface. 


where 


e = q/qb 

q = emissive power of the non-black surface 

q b — emissive power of the black surface 
measured at the same temperature 


(A.1 2) 


The emissivity of the liquid crystal surface is measured using the setup 
shown in Fig. A-l. The blackbody is an aluminum plate with Fresnel rings 
machined into the surface and covered with a black coating. The energy 
radiated from this surface is within 1% of that predicted by Eq. (A.1.1). The 
"radiation thermopile" consists of thirty thermocouples connected in series. 
The hot junctions of the thermocouples are coated with carbon black to absorb 
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the radiation falling upon them and are positioned at the focal point of a 
polished, parabolic reflector. The reference junctions of the "radiation 
thermopiles" are shielded from the incoming radiation and remain at room 
temperature. The EMF produced by the thermopile depends on the difference 
in temperature between the hot and reference junctions, and is a linear 
function of the emissive power of the surface radiating into the thermopile: 

EMF = A + (B)(q) (A.1.3) 

The intercept. A, depends on the surface, while the slope, B, should be 
constant as long as the distance between the surface and the thermopile 
remains the same. 

- The emissivity of the liquid crystal surface is measured by making 
simultaneous measurements of the EMF produced by the thermopile and the 
surface temperature as the surface is heated from room temperature to 100 °C. 
Radiation from the blackbody is used to calibrate the thermopile. For a 
blackbody. 


EMF b = A b + (B)(oT 4 ) (A. 1.4) 

The blackbody is then replaced by a liquid crystal surface. For the liquid crystal 
surface, 

EMFjc = Aic +(B)(e 1c oT 4 ) (A. 1.5) 

or EMFic = Aic + (Be 1c ) (oT 4 ) 


157 



The derivative of EMF b with respect to oT 4 is B, while the derivative of 
EMFic with respect to oT 4 is Be lc . The emissivity of the liquid crystal surface 
can thus be determined by generating curves of EMF vs. oT 4 for both the black 
and liquid crystal surfaces, and taking the ratio of their slopes. 

The value of ej c calculated by averaging the results from the data sets 

was found to be 0.846. The uncertainty in the above is estimated to be 5%. 
Independent measurements of the emissivity made using an emissometer (a 
device for measuring emissivity) yielded readings of 0.85-0.90. 
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A2). Measurement of the thermal cond u ctivity of the lexan/liquid cryst al 
Composite 

The thermal conductivity of the lexan/liquid crystal composite is 
measured using the composite wall shown on Fig. A-2. The wall is 
constructed in a symmetrical manner about a resistance heater patch similar 
to that used by Wang (1984) as a source of constant heat flux. The copper 
plates have grooves cut into them, allowing thermocouples to be placed on 
either side of the composite. The composite is then sealed with epoxy around 

the edges and placed in a water bath. 

The thermal conductivity of the composite is determined by measuring 
the power to the patch heater and the temperature difference across the 
composite. An energy balance on the composite yields 


qi + q2 = qr 


(A.2.1) 


where qj and q 2 represent the heat leaving from either side of the composite 
and q*r represents the power supplied to the heater. An expression for the 
conductivity may be obtained by substituting qi = kATi/Axi, q? = kAT 2 /Ax 2 
into Eq. (A.2.1) and solving for the conductivity, k. The conductivity may 
thus be found if the power supplied to the patch (qr)/ the temperature 
differences across the composite and their thicknesses are known. The results 
of the measurements yielded a thermal conductivity of k=0.1495 W/m/C. 
This compares well with the manufacturers value of 0.146 W/m/C for lexan 
alone. The difference is probably due to the addition of the liquid crystal 
surface and the adhesive transfer tape used to hold the assembly together. 
The uncertainty of the measurement is estimated to be 5%. 
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Figure A.2b-Experimental setup for conductivity measurement. 
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A3). Measurement of a utocorrelation* 

The autocorrelation gives information on the scales and evolution of a 
turbulent flow. It is similar to a frequency spectrum except that the 
information is presented in the time domain rather than in the frequency 

domain. 

The autocorrelation is the correlation of the fluctuating velocity 
component, u’ at two different times, u'(t)u’(t+x). In a steady flow this is 
independent of time, t, but depends only on the tune difference, t. Also, m a 
steady flow the turbulent normal stresses are independent of 
time(u , 2(t)=u’2(t+T)=const.) This can be used to non-dimensionalize the 

autocorrelation as follows: 

. u’(t)u’(t + x) 

pw ~ u' 2 (t) (A.3.1) 

This is known as the autocorrelation coefficient. 

The autocorrelation is related to two important turbulence scales. The 
first is the integral scale, I. This scale is defined as the area under the 
autocorrelation coefficient curve. This time scale represents an average time 
over which u’ correlates with itself. This scale is representative of the large 
scales in the turbulent flow. The second important scale is the Taylor 
microscale, X. This scale is defined by the curvature of the autocorrelation 
coefficient at the origin, as follows: 


•This section was originally written by Mr. Steve Russ. Small modifications 
have been made. 
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2 _ 

X 2 

*=0 K (A.3 2) 

This time scale represents the dissipating scales of the flow. Utilizing the fact 
that the turbulence is stationary, the following relationship can be derived 
(Tennekes and Lumley— 1972): 

du’ ^ 2 _ 2u' 2 
dt ) = .2 

* (A.3.3) 

Thus, the Taylor microscale can be used to estimate the turbulent dissipation 

(assuming small-scale isotropy). Both of these time scales can be converted to 
length scales by multiplying by the local convective velocity, U. 

A simple set of data acquisition and reduction programs to process 
these measurements have been written by Mr. Steve Russ and the author. 
This set-up utilizes a hot-wire anemometer, a Norland Prowler digital 
oscilloscope and an HP lab computer. In this measurement the Norland is set 
to acquire several traces of data from the hot-wire anemometer bridge at fixed 
intervals. The velocity traces are stored on a disk for later data reduction. 
The data reduction program computes the autocorrelation coefficient 
function and the various time scales from the velocity traces. One set of data 
traces with a small acquisition time is required for the measurement of the 
Taylor microscale. The data acquisition rate must be rapid enough so that the 
curvature at the origin is apparent on the autocorrelation curve. A second set 
of data traces with a longer acquisition time is required for the measurement 
of the integral scale. The data acquisition rate must be slow enough so that 
the autocorrelation coefficient curve goes to zero for large t. The 
measurement was accomplished by the following steps: 
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1. A normal hot-wire is placed in the flow and the anemometer is set 
to RUN. 

2. The output from the anemometer is sent to the Norland. 

3. The Norland is set to acquire data at the desired rate. 

4. The high-pass filter is set to 1/2 the acquisition frequency to avoid 

aliasing. 

5. The HP program "DATATRANS" is run. This will take the desired 
data traces. 

6. The HP program "SCALRED" is run to compute the autocorrelation 
coefficient curve and the time and length scales. 

7. The HP program "PLOTRHO" is used to graph the data. 

It was found that at least eight velocity traces at both the high frequency and 
the low frequency were needed to obtain a smooth autocorrelation coefficient 
curve. The rates of data acquisition depend on the particular flow. A listing 
of the programs is attached. 
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221 
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260 

290 
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330 
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360 

370 

380 

390 

400 

410 

420 
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440 
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460 

470 

480 

480 

500 

510 

520 

530 


I PROWLER-COMPUTER INTERFACING PR06RAM (OATATRANS) 

I 

1 TRANSFER OF OATA <-• 

DIM Att 100003 BUFFER ,CS(300) ,A1 (2 ) ,A2<2 ) ,V(2 > ,Uel(2 > 

OIM Fact or (3), Offset <3 > ,Uolt (4095 > ,Uelc< 4095 ) 

REMOTE 715 
! 

I HOT-WIRE INFORMATION 60ES HERE 
I 

A— 1.60786 

8-3.28453 

Powerhw-.435 

INPUT 'INPUT TEMPERATURE OF FLOW* .Temp 
Sqrcf-SQRt 225/ < 250-Temp ) ) 

INPUT 'INPUT ATMOSPHERIC PRESSUREUn. Hfl)'. Press 
Press-Press«25.4 

INPUT 'INPUT 6AIN FROM SI6NAL CONDITIONER' ,6ainhw 
INPUT 'INPUT OFFSET FROM SI6NAL CONDITIONER* , Of fsethw 
INPUT 'INPUT OFFSET OF CHANNEL A (UNIT 71S)*,0ffch1 
INPUT "INPUT BASE FILE NAME* ,8f ileS 
INPUT ’INPUT NUMBER OF DATA SETS* .Nmax 
PRINT ‘DATA SET:* 

! 

FOR Jk-1 TO Nnex 

PRINT Jk " 

FileS»Bfile*&VAL$<Jk ) 

! CREATE DATA FILES 
ASSI6N ©Prowler TO 715 
ASSIGN ©Buf fa TO BUFFER A$ 

I 

I INITIATE DATA TRANSFER 

OUTPUT 715 USIN6 ** ,K* i *_KCGA* ' I SEND CONTENTS OF CH. A IN XFAST BIN 
WAIT .1 

TRANSFER ©Prowler TO BBuf fa (COUNT 8452 
MASS STORAGE IS * :CS80 .700 ,0* 

WAIT 1. 

OUTPUT 715 USING *«,K*»*R* 


! CONUERSION OF BINARY DATA TO DECIMAL VALUES 
! - 
! RESET BUFFER POINTERS: 

CONTROL 0Buf fa, 5} 1 
WAIT 1 . 

! 

! FIND FACTOR AND OFFSET: 

FOR J-1 TO 2 . 

ENTER 6Buf fa USING *t.2A*iCS ! ENTER EXPONENT 

V( 1 )-IVAL(CS.16)< CONVERT FROM HEXADECIMAL TO DECIMAL 

Exp-2*(V( 1 »— 1 28 » ! COMPUTE EXPONENT 

Su-.5 

Power-4. 

FOR K-1 TO 6 

ENTER ©Buffa USING *t,A*(CS ! ENTER FRACTIONAL VALUE 
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540 

V< 1 )«IVAL<C$ . 1 6 )! CONVERT TO DECIMAL 


550 

Tot*V< 1 >/2*Power 


580 

Power*Power+4 . 


570 

5u«Su+Tot 


580 

NEXT K 


590 

IF J-1 THEN ! COMPUTE FACTOR 


600 

Fac t or( 1 >*Su*Exp 


810 

! PRINT 'FACTOR” iNc; *- *»Factor<1 ) 


E20 

END IF 


630 

IF J-2 THEN ! COMPUTE EXPONENT 


640 

Of fset< 1 >«Su«Exp 


E50 

(PRINT "OFFSET” jNc» "* "»0ffset(1 ) 


660 

PRINT 


670 

END IF 


680 

NEXT J 


690 

l ‘ ' 


700 

! ENTER AND I6N0RE REST OF HEADER: 


710 

ENTER SBuffa US1N6 ”*,240A”»CS 


720 

! 


730 

! CONVERSION OF' DATA: 


740 

FOR J«*1 TO 4096 

- 

750 

ENTER. SBuffa USING ”* ,B" »A1 ( 1 )! ENTER' ONE BYTE 


760 

ENTER »Buffa USING "*,B”jA2<T> 


770 

V< 1 )»A2( 1 >*256 .+A1 < 1 > ! TRANSPOSE ORDER OF BYTES 


780 

V< 1 )**V< 1 >-32768 ! OFFSET BY 8000 HEX 

♦ 

790 

V< t )-V( 1 )*Factor( 1 HOffsetf 1 > ! CALCULATE ACTUAL VOLTAGE 


800 

I-J-1 


810 

! PRINT *1 ,V( ” :Nc: " >- " 1 1 ; V( Nc >-0FFCH1 


820 

Volt( I >“V< 1 >— Of f chi 


830 

! 


840 

! LINEARIZE SIGNAL 


850 

! 


860 

Vol t ( I )«< Volt < 1 > /6a i’nhw+Of f set hu >*Sqrcf 


870 

Velc( I >*( A+B* Vol t < I >”2 >"< l/Powerhw> 


880 

IF INT( 1/100 >-1/100 THEN DISP I ,V< 1 >-Offchl ,<>elc<I > 


890 . 

NEXT J 


900 

BEEP 


910 



920 

I 


930 

! STORE DATA IN DATA FILE 


940 

MASS STORAGE IS " : CS80 ,700 , 1 " 


950 

CREATE BOAT FileS,130 


960 

ASSIGN ©Path TO Files 


970 

OUTPUT ePeth;Velc(* > 


980 

MASS STORAGE IS ” :CS80 ,700 ,0“ 


990 

NEXT Jk 


1000 



1 010 

1 


1020 

END 


ORIGINAL PAGE IS 
OF POOR QUALITY 

166 



10 



20 

• THIS PROGRAM IS USED TO CALCULATE THE AUT OCORROL AT I ON FROM 

30 

! VELOCITY MATRICES TAKEN AT FREQUENCIES 2.000 AND 200 H; 

(SCALRED) 

40 



50 

OPTION BASE 1 


E0 

) 

70 

DIM Veil (4095 ) .Rhol (6 .2 ) ,Rho2( 101,2) .Rholeve<6 .2 ) ,Rho2ave( 101 .2 > .Velh( 4095 

INPUT "INPUT NUMBER OF TRACES " ,Nr 


60 

INPUT ‘INPUT H16H/L0W FREQ. FILE BASE NAME' .FilehS .Filel* 


90 

INPUT ‘INPUT HIGH/LOW FREQUENCIES’ .Freqh .Freql 


100 

1 


110 

Rhol ave< 6 ,2 )“0 


120 

FOR K5«1 TO Nr 


130 

Numb$-VALS(K5> 


140 

Nameh$»Fi lehS&NunbS 


150 

Nanel$“Fi lel$&Numb$ 


160 

N* 1 


170 

M*2 


180 

Rhol ( 1 .1 >-l .0 


190 

Rho2( 1 . 1 >-1.0 


200 

Rhol ( 1 ,2 )«*0. 


210 

Rho2< 1 ,2 )«0. 


220. 

j 


230 

! GET HIGH FREQ. VELOCITY MATRIX FROM DISK 


240 



250 

Th*1/Freqh 


260 

Tl*1 /Freql 


270 

MASS STORAGE IS ”:C580 ,700,1" 


280 

ASSIGN 8Path TO NenehS 


290 

ENTER ©Path ; Vel h( * ) 


500 

MASS STORAGE IS " :CS80 ,700 ,0” 


310 

! 


320- 

1 CALCULATE Ubar AND Urns FROM DATA 


350 

I 


340 

Suml «0 . 


550 

Sun2**0 . 


360 

FOR 1-1 TO 4095 


370 

Sunl«Sun1+Velh( I ) 


380 

Sum2-Sun2+Velh' I >*2 


390 

NEXT I 


400 

Ubar 1 - Sun 1 /.4095 


410 

Urns 1 *Sun2/4094-5un 1 *2/4095/4094 


420 

i ' 

- 

430 

! CALCULATE CORRELATIONS AT Tau«1 *T ,Tau-2* t ,Tau«3«T t TAU' 

«4»T ,TAU=S*T 

440 

! 


450 

Sunl “0 . 


460 

Sun2 c 0. 


470 

SunS-O . 


460 

Sur.4 e 0 . 


490 

5un5-0. 


500 

FOR I«1 TO 4090 


510 

IF I / 1 00=* 1NT (1/100) THEN DISP I 


520 

Sunl 'Sunl + ( Ubar 1 -Ue 1 h< I ) )• ( Ubar 1-Vel h( 1 + 1 ) > 


550 

Sun2-Sun2+( Ubar 1 -Vel h( I ) >* < Ubar 1 -Ue 1 h< 1 + 2 ) ) 


540 

Sun3-Sum3+ ( Ubar 1 -Ue 1 h( I ) )* ( Ubar 1 -Ue 1 h< I + 3 ) ) 
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550 SuM-Sun4+( Ubarl -Velh( I > )«(Uber1-Velh( 1+4 > ) 

560 Sum5«Suei5+<Ubar1-Velh< I ) >•( Ubar1-Velh< I+S ) ) 

570 NEXT I 

580 Rhol (2,1 >-Sun1 /4090/Urns 1 

590 Rhol <3,1 >-Sum2/4090/Urms1 

600 Rhol <4 ,1 >-Sun3/4090/Urms1 

610 Rhol (5,1 )-Sun4/4090/Urr>s1 

620 Rhol <5,1 )»Sun5/4090/Urns1 

630 FOR 1-1 TO 5 

640 Rhol (1+1 ,2 )*I«Th 

650 NEXT I 

660 ! 

670 Rho1ave< 1 .1 >-1 . 

680 Rho1ave< 1 .2 >«0. 

690 FOR 1-1 TO 5 

700 Rho1ave< 1+1 ,1 >-Rho1ave( 1 + 1 ,1 )+Rho1 ( 1 + 1 ,1 > 

710 NEXT 1 

720 ! 

730 ! 

740 ! 6ET LONER FREQ. VELOCITY MATRIX FROM DISK 

750 ! 

760 MASS STORAGE IS “ : CS80 ,700 , 1 “ ' , .. 

770 ASS16N «Path TO Namel$ 

780 ENTER §Pathi Vel 1 < • ) 

790 MASS STORAGE IS “ :CS80 ,700 .0" 

800 ! 

810 ! CALCULATE Urms AND Ubar 

720 ' 

830 Suwl-0. 

640 Sun2*0. 

850 FOR 1-1 TO 4095 

860 SuM-Suml+Vel 1 < I > 

870 Sun2-Sur,2+Vel 1 < I ) "2 

880 NEXT I 

890 Ubar2*Sur r >1 /4095 

900 Urn>s2=Su*-.2/4094-Sum1 '2/4095/4094 

910 ! 

920 ! CALCULATE CORRELATIONS 

930 ! 

940 FOR 1=1 TO 100 

950 Sur.1-0. 

960 FOR J-1 TO 3995 

970 Sunl-Sunl +< Ubar2-Vel 1 < J > )•< Ubar2-Vel 1 < J+I > ) 

980 NEXT J 

990 Rho2< 1+1 ,1 >«Sum1 /3995/Ur-ns2 

1000 Rho2< 1 + 1 ,2 )«I«T1 

1010 DISP I 

1020 NEXT I 

1030 Rho2ave< 1 ,1 >“1 . 

1040 Rho2ave< 1 ,2 )-0 . 

1050 FOR 1-1 TO 100 - 

1060 Rho2ave< 1 + 1 ,1 >-Rho2ave< 1 + 1 ,1 >+Rbo2< 1 + 1 ,1 ) 

1070 . NEXT 1 

1090 NEXT K5 
1050 FOR 1=1 TO 5 
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1100 Rho1ave(I+1,1)“Rho1ave(I+1.1)/Nr 

1110 Rholavef 1 + 1 .2 >*Rho1 ( 1 + 1 ,2 ) 

1 120 NEXT I 

1130 FOR 1 = 1 TO 100 

1140 Rho2ave< 1+1,1 >=Rho2eve< 1+1 ,1 >/Nr 

1150 Rho2ave( 1 + 1 .2 >«Rho2< 1 + 1 .2 ) 

1160 NEXT I 

1170 ! *.*****••*••«••*••*•*•«••••••••* 

1180 ! 

1190 ! LEAST SQUARES FIT PARABOLA TO THIS DATA TO GET MICROSCALE 

1200 ! 

1210 Sun 1 «0 . 

1220 Sun2=0. 

1230 Sun3=0. 

1240 Sun4«0. 

1250 FOR 1 = 1 TO 6 
1260 Y-Rho1ave< I ,1 > 

1270 X“Rho1ave( I ,2 >“2 

1280 Sun1=Sun1+X 

129G Sun2“5un2+Y 

1300 Sun3“Sun3+X«X 

1310 Sun4=Sun4+X» Y 

1320 NEXT I 

1330 Slope=( 6*Sun4-Sun1 «Sun2 >/(6«Sun3-Sun1 *Sun1 ) 

1340 Microt=SQR(-1 .0/Slope ) 

1350 Microl=Ubar 1 +Microt 
1360 ! 

1370 ! CALCULATE INTEGRAL SCALES 
1380 ! 

1390 Sun 1 =0 . 

1400 FOR 1=2 TO 100 STEP 2 

1410 Sunl-Sunl+(Rho2eve< 1-1,1 1+4 . 0«Rho2eve< I ,1 )+Rhc2eve( 1 + 1 ,1 ) )•( RhoZeve! 1+1 , 

2 >-Rhc2ave< 1-1 ,2) )/6 

1420 NEXT I 

1430 Itine=Sun1 

1440 I lengt h=Uber2«It ine- 

1450 ! 

1460 i PRINT OUT RESULTS 

1470 INPUT 'DO YOU WISH TO STORE DATA? <Y/N)",StoS 
1480 IF StoS'" Y" THEN 

1490 INPUT "INPUT DATA FILE NAME" ,Stor$ 

1491 MASS STORAGE IS ":CS80. 700.1" 

1500 CREATE BOAT StorS.7 

1510 ASSIGN 8Peth TO Stor$ 

1520 OUTPUT 0Path;Rho1ave( • > ,Rho2ave< * ) 

1521 MASS STORAGE IS " : CS80 ,700 ,0“ 

1530 END IF 

1540 ! 

1550 INPUT "DO YOU WISH RESULTS PRINTED ON SCREEN OR PRINTER ? <S/P>'.Pri* 

1560 PRINTER IS 1 

1570 IF Pr j $= ” P " THEN PRINTER IS 701 

1580 T i2=SQR< Urn 5 2 )/Ubar2« 100 

1590 PRINT- "AVERAGE VELOCITY " ,Uber2 

1600 PRINT "AVERAGE TURBULENCE INTENSITY ”.Ti2 

1610 PRINT "INTEGRAL TIME SCALEis) ’.IUne 
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1620 PRINT " INTE6RAL LEN6TH SCALE <m> ’.I length 

1E30 PRINT "TAYLOR MICRO TIME SCALE (s) " .Microt 

1 640 PRINT "TAYLOR MICRO LENGTH SCALE <m> " ,Micr-ol 

1650 PRINT " TAU RHO " 

1660 FOR 1-1 TO 11 

1670 PRINT Rho2( ( I — 1 )»10+1 ,2 ) ,Rho2< (1-1 )*10+1 ,1 ) 

1680 NEXT I 
1690 ! 

1700 PRINTER IS 1 

1710 ! 

1711 ! CALCULATE LEN6TH SCALES FROM TIME SCALES 

1712 ! 

1713 INPUT "DO YOU WISH TO CHANGE TO LENGTH SCALES?" ,Inp$ 

1714 IF Inp$=" Y" THEN 

1715 FOR 1-1 TO 101 

1716 Rho2< 1 ,2 )-Rho2( 1 ,2 >«Ubar2 

1717 NEXT I 

1718 FOR 1=1 TO 6 

1719 Rhol ( I .2 )-Rho1 ( I ,2 )«Uber1 

1720 NEXT I 

1 722 ! 

1723 INPUT "DO YOU WISH TO STORE DATA? <Y/N)",StoS 
1728 IF 5to$=“Y" THEN 

1738 INPUT "INPUT DATA FILE NAME" .StorS 

1739 MASS STORAGE IS " : CS80 ,700 , 1” 

1748 CREATE BDAT StorS ,7 

1758 ASSIGN SPeth TO StorS 

1768 OUTPUT ®Path;Rho1ave( • ) ,Rho2eve( * ) 

1769 MASS STORAGE IS 1 :CS80 ,700 .0" 

1778 END IF 

1788 ! 

1798 END IF 
1808 i 
16)8 END 
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1C I ***** 

20 ! THIS PRQ6RAM IS USED TO PLOT DATA EITHER MANUALLY OR THR0U6H A DATA FILE 

30 ! ON LINEAR-LINEAR AXES ( PLOTRHO ) 

40 * ♦ 

50 OPTION BASE 1 

60 DIM Xd<500) .Yd< 500). TitleSt 50 3 ,Labelx$[50) ,Labely$C50 3 ,Rho2( 100,2 ) ,Rhc1<£, 

2)! ARRAY TO BE PLOTTED 

70 GRAPHICS ON 

80 6CLEAR 

90 6INIT 

100 LOR6 5 

110 DE6 

120 * * 

130 Title*-* 4 " ! TITLE OF PLOT 

140 Xmin«0. ! MINIMUM UALUE OF X 

150 Xroax* . 02 ! MAXIMUM UALUE OF X 

150 Yni n»- . 2 ! MINIMUM UALUE OF Y 

170 Ynax-1. ! MAXIMUM UALUE OF Y 

180 Xt ic* . 002 • SMALL SCALE 

190 Nxt ic “2 ! HOW MANY SMALL SCALES IN LARGE SCALE 

200 Yt ic* . 05 

210 Nytic*4 

220 Leb$* " Y" ! WANT LABEL 

230 LabelxS= M T [S3" ! X-AXIS LABEL 

240 Label y$* "RHO" ! Y-AXIS LABEL 

310 * ♦ i 

320 LDIR 0 . 

330 CSIZE 5 

340 LORG 5 

350 FOR I*-. 1 TO .3 STEP . 1 

350 MOUE 70+1,95 

370 LABEL Titles 

390 NEXT I 

390 CSIZE 5 

400 LORG 5 

410 MOUE 59,5 

420 LABEL Label*® 

430 LDIR 90 

440 MOUE 5,52 

450 LABEL LabelyS 

450 UIEWPORT 15,124,12,90 

470 ! 

480 ! LINEAR-LINEAR AXES 

490 WINDOW Xmi n ,Xmax , Yru n , Ymax 

500 AXES Xt ic ,Yt ic ,Xmin f Ymin f Nxt ic ,Nyt ic ,5 

510. AXES Xt ic , Yt ic ,Xnax .Yna *. ,Nxt ic ,Nyt ic ,5 

520 • IF Grid®* “Y’* THEN 6RID Nxt ic*Xt ic ,Nyt ic* Yt ic ,Xmax , Ymax 

530 ! 

540 CLIP OFF 

550 LDIR 0 

550 IF LabS* "N" THEN 60T0 800 

570 Ylab*( Yma*-Ymin )/20 

580 FOR 1*1 TO < Xma * - Xri i n ) / < X 1 1 c • N>; t i c ) 

590 X 1 * I «Xt ic*N>:t lc+Xmi n 

600 MOUE Xl.Ynin 
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610 L0R6 6 

620 LABEL XI 

630 NEXT I 

640 I 

650 Xlab=( Xroax-Xrnn )/25 
660 LOIR 90 

670 FOR I-T TO ( Yrnax-Ymin )/( Ytic«Nytic >+l 
680 Yl-I«Ytic«Nytic+Ymin . 

690 HOWE -Xlab+Xnin ,Y1 

700 IF AB5(Y1 XI .E-10 THEN GOTO 720 

710 LABEL Y 1 

720 IF ABS(Y1 XI .E-10 THEN LABEL '0" 

730 NEXT I 

740 LDIR 0 

750 CLIP ON 

760 MOVE 0,0 

770 LINE TYPE 5 

780 IDRAUI 100 ,0 

781 LINE TYPE 1 

790 * 

800 ! PLOT DATA ENTERED MANUALLY 

810 CLIP ON 

820 LORG 5 

830 CSIZE 1 , . 5 

840 INPUT "INPUT DATA FILE NAME" .FileS 
850 MA55 STORAGE IS ":CS80, 700.1" 

850 A5SI6N ©Path TO FileS 

870 ENTER 0Path;Rho1<« ),Rho2<« > 

880 MASS STORAGE IS * :CS80 ,700 .0" 

890 FOR 1=1 TO 100 

900 MOVE Rho2< 1 ,2 >,Rbo2< 1 ,1 > 

910 LABEL “+" 

920 NEXT I 

921 FOR 1=1 TO 6 

922 MOVE RhoK 1 ,2 I.RhoK T,1 > 

923 LABEL "+" 

924 NEXT I 

930 * * 

940 INPUT 'DO YOU WISH A HARD COPY? (Y OR N)”,A5$ 

950 IF A5$=“N" THEN 1040 

960 INPUT "EXPt.r'OED MODE? < Y/N )" .Expands 

970 IF Expands*” Y" THEN 

960 DUMP DEVICE IS 701 .EXPANDED 

990 DUMP 6RAPH1CS 

1000 END IF 

1010 IF E*pandS*‘N" THEN DUMP 6RAPHICS 1 TO *701 
1020 i OUTPUT 701 ; " 

m 

1030 OUTPUT 701 i" 

1040 GCLEAR 
1050 END 
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